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HURRICANE EFFECT ON SEA LEVEL AT CHARLESTON 


Bernard D. Zetler* 
(Proc. Paper 1330) 


ABSTRACT 


Tropical cyclones from 1922 to 1955 have been classified according to the 
orientation of their tracks with respect to Charleston Harbor. Theoretically, 
the displacement of sea level depends on this orientation; the tabulated dis- 
placements caused by tropical cyclones during 34 years of tide observations 
bear out the relationships. 


When it was decided in the Coast and Geodetic Survey to make a study of 
the effect of hurricanes on sea level, Charleston was chosen as the area of 
the pilot study for several reasons. Its position makes it vulnerable to peri- 
odic hurricanes, and the low elevation of its downtown section creates the 
need for accurate forecasting of storm heights. Furthermore the problem 
appeared to be simpler for a port just off the ocean than, for example, a port 
on Chesapeake Bay in which interference oscillations may be anticipated. 
Lastly there was available for Charleston a relatively long series of tide 
observations dating back to 1922. 

The attempt in this study has been to classify hurricanes into categories 
such that some groups are likely to have a large effect on sea level while 
others have little or none. Figure 1 shows the basis of one set of criteria 
which was chosen. Considering the Atlantic coastline near Charleston to lie 
in a SW to NE direction, the optimum direction of water motion to raise sea 
level is from the SE. Considering the Coriolis effect in the northern hemi- 
sphere of water moving about 45 degrees to the right of the wind, water move- 
ment from the SE will be set up by easterly winds. If an observer faces into 
the wind, the storm center will be on a line about ten points (112°) to the 
right. Therefore with easterly winds at Charleston, the storm center will 
have a bearing of approximately SSW. Accordingly, all other things being 


Note: Discussion open until January 1, 1958. Paper 1330 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 1, August, 1957. 


*Oceanographer, Coast and Geodetic Survey, U. S. Dept. of Commerce, 
Washington, D. C. 
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equal, storms having a bearing from Charleston of about SSW will cause a 
rise in sea level at Charleston while storms having a bearing of about WNW 
or ESE will cause currents to move parallel to the coast with little effect on 
sea level. A tropical cyclone bearing NNE from Charleston should lower 
sea level by moving water away from the coast. Storms having other bear- 
ings from Charleston will have effects proportionate to their departures 
from the above directions. 

Figure 2 shows eight octants radiating out from Charleston, each repre- 
senting an arc used in classifying the bearing of a cyclone center. Each 
name represents the direction in the center of the octant. By tabulating the 
octants through which the center of a cyclone moves in approaching and mov- 
ing away from Charleston, one really is enumerating the changing bearing of 
the storm center as observed from Charleston. 

Those storms classified as having Charleston on the right side of their 
tracks ordinarily move inland whereas those with Charleston on the left side 
remain out at sea. The criteria was first considered important on the as- 
sumption that storms moving inland pass a critical depth in shallow water at 
which resonance contributes to a surge.(1) Using this reasoning in studying 
the surge effect at Charleston, it would be necessary to distinguish between 
storms crossing the Atlantic coast and those crossing the Gulf coast. How- 
ever this study shows large positive displacements with the track crossing 
the coastline in the Gulf of Mexico as well as the Atlantic coast. Apparently 
the criteria of Charleston being on the right side of the track is important 
because of associated strong winds over shallow water outside the harbor 
rather than because of a surge being set up at a critical depth. A storm with 
Charleston on the left side of its track contributes to changing sea level in 
the form of forerunners which are moderate rises in sea level over a period 
of time due largely to currents set up by wind circulation. In deep water 
much of this effect is dissipated before reaching the coast. 

The question of how close a cyclone had to come to Charleston to affect 
sea level appreciably was resolved by including all storms shown on the 
track charts for the first few years considered. It became apparent that it 
was possible to restrict the study to cyclones coming within 345 miles (300 
nautical miles) of Charleston. This is a convenient figure to work with as it 
is measured on a chart by five degrees of latitude near the latitude of 
Charleston. In the final tabulations, only the portions of tropical cyclones 
within a radium of 345 miles of Charleston are considered. This probably 
includes some cyclones which need not be considered but the effective 
diameters of tropical cyclones vary considerably in size and the diameters 
of particular storms are not usually listed in the published statistics. 

The cyclones are classified first according to the octants in which they 
came closest to Charleston but the data on other octants are also important. 
For example two cyclones may have been closest to Charleston in Octant ESE. 
If one of them has passed through Octants SSW and SSE in approaching ESE, 
it presumably is more likely to have raised sea level than the other. The 
final tabulation also subdivides the storms classified within an octant to 
those with Charleston on the right side of the cyclone track and those with 
Charleston on the left. Tropical cyclones with Charleston on the right side 
are presumed to have a greater effect on sea level. The results bring this 
out very clearly. 
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WINDS 


OPTIMUM WATER 
MOVEMENT 


Figure 1. Theoretical bearing of center of tropical cyclone for maximum 
rise in sea level at Charleston. 


Figure 2. Octants radiating from Charleston used in classifying 
the bearing of the center of a tropical cyclone. 
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Tide Records 


In determining the displacement of water level caused by a tropical cy- 
clone, it was found desirable and necessary to consider the continuous record 
as represented by the complete tide curve. In most cases, a comparison of 
the observed and predicted heights on each hour presented a true picture, but 
it was found unsatisfactory for final results to compare heights only at the 
times of high and low tides, for, being so far apart in time, they seldom 
showed the maximum variation. Furthermore the time of high or low tide 
may be advanced or retarded several hours by a storm and a comparison of 
observed and predicted high or low tides would deal with water level at dif- 
ferent times. 

Despite these limitations comparisons of observed and predicted high and 
low tides were used as the preliminary step in determining what storms were 
to be included in the study. Any storms for which the calculated displace- 
ments did not include at least one value of one foot or more were eliminated 
from further consideration. These are shown in the tabular results as having 
a displacement of less than one foot. Obviously, in using this arbitrary rejec- 
tion limit, some storms in which the calculated maximum hourly displace- 
ment would have been somewhat over a foot were rejected from further study. 
However, in each case the observed tide curve was examined for unusual 
fluctuations between the high and low waters in order that no important dis- 
placements would be rejected by this method. 

Predicted tides at Charleston are referred to a datum of mean low water, 
based on gage records for the 19-year period 1924 - 1942, and take into ac- 
count seasonal variations in sea level. Inasmuch as sea level along our east 
coast has been rising gradually (approximately 0.4 foot in the last few dec- 
ades), each year’s observations were referred to mean low water observed 
in the same year before being compared with predictions. 


Tabulation of Tropical Cyclone Data and Sea Level Displacement 


The following tables furnish data on tropical cyclones from 1922 to 1955, 
describing their paths, speeds, and their bearings and distances from 
Charleston. They indicate the maximum displacement of sea level and the 
interval of time between the cyclone coming closest to Charleston and the 
maximum displacement. When the maximum displacement occurred before 
the cyclone reached its closest point, the time lag is preceded by a negative 
sign. The maximum water height (referred to Mean Low Water) attained 
during the period of the storm is given in the last column. Underlined names 
of storms indicate that they were hurricanes. 

Tropical cyclones are numbered in accordance with Weather Bureau track 
charts except that those earlier than 1933 are numbered in accordance with 
the track charts in “HURRICANES” by I. R. Tannehill. Numerous adjectives 
were found to describe the intensity of tropical cyclones in various sources. 
Even the clear-cut distinction of whether or not a tropical cyclone attained 
hurricane status is somewhat unsatisfactory since that degree of intensity 
may not apply to the storm during the period it may have affected sea level 
at Charleston. Some track charts show the position of the storm center 
every twelve hours; others give it every twenty-four hours. It was necessary 
to infer that the center of the storm moved at a uniform rate between adjacent 
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TABLE 1 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant SSE When Closest to Charlestm 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 
Cyc lone Octant Fath Duration Speed Charleston Height La MLW 
miles hours m.peh. miles feet hours feet 


Charleston on right side of cyclone track 


1940 II SSE 110 11 10 1,0 +2,2 
SSW 235 18 13 170 


1951 Able ESE 215 iu 20 235 +1,1 
Approach) SSE 220 20 230 
Charleston on left side of cyclone track 
1938 VII WSW 4 51 
SSW 50 
(50) 
(50) 
(49) 


10 
16 
13 
28 
25 


26 


3 
(0.4) 
(0.4) 
(7) 


7 
19 
9 
3 
2 
iu 
u 
3 
4 
le 


1953 Hazel 


=) 
VAD 


1938 VI 
1941 III 
(ist Approach) 


1931 IV 


wF 


»t 


ASCE 1330-5 
-2 59 
- 6.8 
1924 IV WSW 185 175 <1,0 - 66 
SSW 145 55 
SSE 40 45 
ESE 55 50 
ENE 275 70 
1932 IV SSW 285 105 <1.0 - 6h 
SSE 65 22 70 
ESE 95 2h 75 
ENE 270 23 115 
1937 1 ssw 240 17 145 <1.0 - 505 
SSE 105 18 85 
ESE 105 1s 90 
ENE 265 15 115 
ssw 110 22 265 41.3 =22 
SSE 195 28 205 
L ESE 215 27 215 
SSE up | ? 230 <1.0 6.9 
ESE L5 15 295 
ssw 10 275 <1.0 62 
SSE 195 9 250 
ESE 100 ? 270 
1948 VIII SSE 195 a 18 320 <1,0 « 2 
ESE 100 25 325 


1330-6 HY 4 August, 1957 


TABLE 2 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant SSW When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 
Cyclone Octant Path Duration Speed Charleston Height La MLW 
miles hours m.p.h, miles feet hours feet 


Charleston on right side of cyclone track 


ESE 


SSE 250 32 35 
SSW 


WNW 


1947 VIII 


1936 Ix 


1937 III 


1939 II 16 13 
(1st Approach) WSW 40 5 8 


Charleston on left side of cyclone track 


1946 VI SSW 285 


1950 XII SSW 14,0 


1937 VI WSW 90 
85 


1941 VI SSW 220 ? 


? 60 +14 2? 
- 230 +167 2 607 


1.0 6.5 


- 2L0 < 1.0 720 


* According to measured high water mark, At time gage ceased operating, 
displacement was +3.7 feet, 


1940 III 13 8 280 +5,1* 2 10.7 
wsw 95 15 6 35 
250 23 11 125 
1947 IX SSE 95 6 16 315 42.6 48S 
ESE ) Course ( 6 ba (300) 
Juncertain(” ye (215) 
SSE 75 8 9 80 
SSW 70 6 12 70 
wsw 95 ? 100 
SSE 210 (13) (16) 160 41.3 6.7 
ssw 150 7 21 up 
WSW 200 ? 160 
SSE 125 6 2 27 <21.0 69 
SSW 12 128 24,0 
WSW 125 ( 9) (14) 275 
SSE 85 a 300 <1,0 6.6 
SSW 230 22 10 250 
wsw 145 yu 10 275 
330 
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TABLE 3 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant WSW When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Length Distance of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height Lag MLW 
miles hours’ m.p.h. miles feet hours feet 
Charleston on right side of cyclone track 


1953 Florence WSW ? ? -72 707 


1934 I =3 


om aw 


3 
(12 
(4 
18 
12 
3 


+2.3 6.9 


1933 VI <1,0 6.6 


1953 Alice ? <1,0 569 


1948 III WSW 5 <1.0 73 
WNW 


* The maximum displacement may have been exceeded and the time lag changed 
during the period the gage was inoperative, 


ay 
ASCE 1330-7 
SSE 145 16 130 
WSW 13 LO 
WNW 45 8 50 
NNW 80 5 115 
WNW 150 - 95 
1927 V SSE 145 18 205 (+1.9 +1 5-3 
ssw 150 15 80 (-1.0 +8 4 
WwSW 50 17 60 
WNW 65 16 65 
NNW 225 ) (19) 85 
NNE 65 ) (16) 270 
1952 Able SSE 205 190 +3,0% +2% 6.6 
SSW 155 13 70 
viSW 55 18 60 
WNW 70 18 65 
NNW 16 13 105 
NNE 80 ( 6) (13) 270 
1950 SSW 60 180 42.7 +4 6.8 
wsw 190 13 15 175 
WNW 85 14 270 
1926 I SSW 255 13 20 215 42.8 +2 707 
wsw 10 2 210 
‘ WNW 60 3 20 265 
1950 XI SSW 18 13 230 
wsw 230 ? 225 
wsw 180 ? - 235 
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Tropical 


1924 V 
WNW 


1947 VII SSW 


1946 V SSW 
WNW 
NNW 
NNE 
ENE 


SSW 
WSW 
WNW 
NNW 
NNE 


Length 
of 


miles 


hours 


The tide gage at Charleston was 


Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant WNW When Closest to Charleston 


Storm Data for Octants Traversed 


Minimum 
Distance 


to 


Charleston on right side of cyclone track 


11 
1 
1 

(10) 


30 
25 
25 
35 


110 


Duration Speed Charleston Height 
feet hours 


(41.9 
(-2.2 


+242 


inoperative during this storm, 


* The maximum displacement may have been exceeded and the time lag changed 
during the yseriod the gage was inoperative. 


Maximum 
Displacement 
Of Sea Level 


August, 1957 


Maximum 
Hei ght 
Above 


Cyclone Octant Fath Pp g z£ MLW 
miles 


feet 


6.7 


6.5 


702 


+1,0* (-40)# (6.3) 


HY 4 
TABLE ha 
320 2 
NNW 25 25 
NNE 310 (31) 
1945 IX SSW 260 26 10 | | || -3 7.9 
WwSW 75 5 15 60 
WNW 50 3 17 50 
NNW 90 6 15 6 
NNE 205 12 17 135 
1944 X SSW 17 15 105 4308 <3 8.5 
WSW 85 4 21 65 
WNW 45 2 22 60 
NNW 100 5 2 65 
NNE 240 11 22 130 
WSW 130 19 75 
WNW 45 22 60 
NNW & 20 65 
NNE 275 21 95 
215 10 22 150 +2,3 =, || 
105 5 21 90 
75 3 25 85 
100 8 12 90 
205 15 4 125 
100 I 6) 17 270 
1935 II 
1933 XII SSW 200 28 7 270 Pp 
WSW 230 2h 10 125 
80 ? 100 
NNW LO ? 105 


TABLE 4b 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant WNW When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Distance Of Sea Level Height 
Tropical to Time Above 
Cyclone Octant Fath Duration Speed Charleston Height Lag MLW 


miles hours my.h. miles feet hours feet 


Charleston on right side of cyclone track - cont'd 


1949 II 185 13 (41.9 +5 765 


165 15 (-1.5 +10 
pate) 22 


230 26 
50 


1928 II WSW 235 12 
WNW 260 12 
NNW 100 
1932 VII WNW 


The tide gage at Charleston was inoperative during this storm, 
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4 ~ 
WSW 
WNW 
NNW 5 
NNE 305 
310 <1,.0 = 6.4 
300 
310 - - - 
NNW 
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TABLE 5 


August, 1957 


Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant NNW When Closest to Charleston 


Storm Data for Octants Traversed 


of 


Minimum 
Distance 
to 


Maximun 
Displacement 


Of Sea Level 


Time 


Maximum 
Height 
Above 


Cyclone Octant Fath Duration Speed Charleston Height Lag MLW 


1939 II 
(2nd Approach) 


1934 II 


miles 


hours 


MePehe 


miles 


feet hours 


Charleston on right side of cyclone track 


SSW 
WSW 
WNW 
NNW 
NNE 


WSW 


295 


19 


6 
2 
2 
2h 
26 
5 
2 
2 
9 
18 
2h 
5 
5 
12 
32 
2h 
12 
10 
10 
6 


o 


16 
15 


+21 


feet 


(8.1) 


HY 4 
i 
= 
1928 V = 90 ( +3.8 
90 30 (-1.2 
20 10 20 
20 10 15 
330 20 
1941 V | 12 70 <1,0 567 
WNW 70 uu 35 
NNW 25 12 20 
NNE 30 15 25 
" ENE 115 13 35 
ESE 235 13 130 
1929 II wsw 310 13 90 (+1.6 +L, 720 
WNW 70 65 (-2.2 +10 
NNW 80 16 & 
NNE 250 21 100 
1928 I SSW 220 7 285 6.3 
WSW 230 10 225 
WNw 155 13 155 
NNW 1,0 uy 145 
NNE 200 20 190 
1936 XVI WSW 170 15 24,0 <1,.0 - 6.3 
WNW 170 12 165 
NNW 160 27 160 
7 NNE 160 27 220 
WNW 200 313 270 <1,0 - 5.7 
NNW 230 19 250 
NNE 35 i 3 310 
NNw 20 310 <1.0 - 503 
#2 
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TABLE 6 


Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant NNE When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 

Length Distance Of Sea Level Height 

Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height Lag MLW 
miles hours mepeh. miles feet hours feet 


Charleston on left side of cyclone track 


1942 VI ENE 50 4 12 310 <1.0 - 720 
NNE 115 240 


points to obtain a time lag measured in hours. Occasionally an extrapolation 
was necessary for this measurement and parentheses are set around these 
figures to indicate an inference. A number of cyclones dissipated near the 
points on their tracks closest to Charleston; a time lag estimate would be un- 
reliable for these and a question mark appears in the column. Within each 
classification, the tropical cyclones are listed in order of increasing mini- 
mum distance to Charleston. Several cyclones came within the 345 mile 
range of Charleston, moved away, and then came back again. These ap- 
proaches are listed independently in the tables but are distinguished by a 
notation “lst approach” or “2nd approach.” 


Tropical Cyclones In Octant SSE 


Theoretically a cyclone in Octant SSE will raise sea level but not as much 
as one in SSW. However a cyclone approaching Charleston from SSW will 
usually pass over Florida and in doing so, lose some of its energy. 

There are 11 storms listed in this classification. There was a rise in 
sea level for the two cyclones with Charleston on the right side of their 
cyclone tracks. All nine having Charleston on the left showed little dis- 
placement. 

Hurricane II, 1940, had a displacement of +2.2 feet 24 hours after coming 
closest to Charleston. This storm moved from Octant SSE into Octant SSW 
and the maximum displacement took place while the storm was in SSW, ac- 
counting for the large time lag. 

Hurricane Able, 1951, had a maximum displacement of +1.1 feet, two hours 
before coming closest to Charleston in its first approach. It then turned 
counter-clockwise, out of range of Charleston, without crossing the coast, 
thus accounting for the small displacement and small lag. Hurricane II, 1940, 
came almost a hundred miles closer to Charleston than did Hurricane Able, 
1951, and moved considerably slower. 

Four of the remaining nine came within less than one hundred miles of 
Charleston and all four passed through Octant SSW. One of them, IV 1924, 
attained hurricane intensity but had a small diameter. The fact that 
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TABLE 7 


Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant ENE When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacemannt Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 
Cyclone Octant Path Duration Speed Charleston Height Lag MLW 


miles hours mepeh. miles feet hours feet 


Charleston on left side of cyclone track 


1954 Hazel SSE 220 8 28 45 a7 
ESE 105 4 26 100 (-1.5 15 
ENE 105 3 35 90 
NNE 255 5 51 125 

1944 ITT SSE 165 15 180 és 
ESE 140 9 16 105 
ENE 8 5 16 100 
NN 230 13 120 


10 165 
32 6 160 
4 15 315 
15 185 +1.7 +5 7.8 
12 12 170 
17 4 185 


8 290 <1,0 - 5.8 
25 9 185 
15 13 180 

3 12 265 


220 


1941 III 
2nd Approach) 


1955 Connie ESE 255 
ENE 200 
NNE 
1955 Diane ESE 215 
ENE 145 
NNE 230 
1953 Barbara SSE 
ESE 230 
ENE 195 
NNE 35 
1955 Ione ESE 29 13 1s = 6.9 
ENE 160 13 12 215 
NNE 90 10 9 255 7 
1944 VI ESE 230 wu 16 230 <1.0 - 6.1 
ENE 235 10 2h 225 
1933 XIII ESE up 15 9 270 <1.0 = 509 
ENE 195 18 lu 265 
1936 XIII ENE 185 16 12 310 18 - 6.1 
1940 IV ESE 35 4 9 330 <1,0 - 6.9 
ENE 145 15 10 310 
1933 VIII ENE 105 8 13 325 41.7 =2 6.6 
NNE 50 4 12 330 . 
ESE 50 6 8 330 4.8 <2 To2 
ENE 100 lu 9 325 
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Charleston was on the left side of their hurricane tracks appears to be the 
explanation for their failure to raise sea level appreciably. 


Tropical Cyclones in Octant SSW 


Theoretically, cyclones with centers in Octant SSW will have the maximum 
ability to raise sea level at Charleston. Of the six cyclones which had 
Charleston on the right side of their tracks, the three closest caused notable 
rises in sea level; the three furthest away did not. Furthermore the amount 
of displacement by the first three appears directly related to their minimum 
distances to Charleston. This ordered sequence may also be due to the first 
two cyclones, III 1940 and [IX 1947, attaining hurricane intensity while the 
third, VIII 1947, did not. The time lags for the second and third cyclones are 
negligible, as is to be expected of most cyclones coming closest to Charleston 
in Octant SSW. The time lag on the first, with a displacement of +5.1 feet, 
was not determined because of tide gage difficulties. 

The loss of some tide data for Hurricane III 1940 illustrates a basic diffi- 
culty in a study such as this. The storms which strike Charleston hard and 
cause large displacements in sea level frequently cause tide gage failure and 
on a few occasions have destroyed the tide gage installation. The loss of 
record of some extreme conditions is a serious handicap to a study such as 
this. The program now being carried out of modifying the tide gage design, 
raising the gages and remoting recorders will help in the future. 

There were four cyclones which were closest to Charleston in Octant SSW 
with Charleston to the left of their tracks. Tropical cyclone VI 1946, which 
came within 60 miles of Charleston caused a displacement of +1.4 feet. Hur- 
ricane XII 1950, was 230 miles from Charleston and caused a displacement of 
+1.7 feet. The other two, both 240 miles from Charleston, had little effect on 
sea level at Charleston. Since all four dissipated in Octant SSW there is no 
data on time lag. 


Tropical Cyclones in Octant WSW 


As is to be expected in view of the usual paths of West Indian tropical cy- 
clones, all the storms which came closest to Charleston in Octant WSW have 
Charleston on the right side of their cyclone tracks. Since many of these 
passed through Octant SSW and even SSE prior to reaching WSW, a large 
positive sea level displacement with little time lag is to be expected from 
these cyclones. This is found to be the case. 

Six of the ten cyclones in this category raised sea level from 2 to 4 feet; 
the time lag in each case was small. Of the remaining four, one dissipated in 
the octant and the other three were not hurricanes and their distances from 
Charleston were large. 

Sometimes the elevation of sea level, as caused by a tropical cyclone, is 
followed by a lowering of sea level appreciably below the predicted value. 
The maximum negative displacement sometimes occurs just a few hours after 
the maximum positive displacement; at other times the time difference is a 
day or more. As noted previously, tropical cyclones in Octant NNE should 
lower sea level at Charleston most effectively and cyclones in Octants NNW 
and ENE should have similar effects. None of the tropical cyclones classi- 
fied as coming closest to Charleston in Octants SSE and SSW passed through 
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TABLE 8a 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant ESE When Closest to Charleston 


Storm Data far Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height Lag MLW 
miles hours mph. miles feet hours feet 


Charleston on right side of cyclone track 


1934 III 120 ? & “1.2 -164 5.6 
50 8 50 
up 16 60 
255 17 


100 10 
230 23 
SSE 70 7 


Charleston on left side of cyclone track 


160 
110 


BERE 


21 
6 
6 


1922 III 


1935 265 41.4 413 bob 
300 
1930 IT ssw <1,0 - 508 
» 
ESE 115 10 100 
ENE 235 13 1,5 
1945 I SSW 13 uo 42.6 +1 7.0 
| SSE 100 17 105 
| ESE 10 17 100 
ENE 265 19 125 
1925 III ssw 16 6 27 230 +2,0 +2 6.9 
SSE 165 n 15 120 
ESE 110 6 18 110 
ENE 165 13 13 130 
NNE 110 9 12 230 
1954 Carol SSE 80 17 5 310 +20 #4 6.8 
ESE 260 72 L up 
ENE 255 15 17 160 
< 1934 VI SSE 190 17 iu 185 <1.0 ~ 509 
ESE 1,0 6 23 150 
ENE 230 8 29 170 
SSW & ? 290 +1,.8 702 
SSE 195 (2) (7) 185 
ESE 190 39 5 180 
ENE 100 18 6 260 
1951 How ssw 50 3 17 320 +7 
SSE 215 12 18 215 
ESE 180 12 15 195 
ENE 150 18 8 270 
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TABLE 8b 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant ESE When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimun Displacement Maximum 
Length Distance Of Sea Level Heigtt 
Tropical of to Time Above 
Cyclone Octant Path Duration Speed Charleston Height La MLW 
miles hours mph. miles feet hours feet 


Charleston on left side of cyclone track - cont'd 


SSE 15 0.7 330 +163 (+14) 69 
ESE 205 ? 200 


SSE 110 13 
ESE 190 18 
ENE 160 12 


<1,0 6.7 


SSE 195 
ESE 195 
ENE 185 


6.1 


ESE 
ENE 


7.0 
6.8 
506 
6.9 


602 


Tide gage not oper= 
ating at this time, 


< 1,0 6.7 


(2nd Approach) 


1943 II <1.0 


1932 IIT <1,0 


BR uF BE &e wo who BB 


21 

8 
u 
13 

5 
16 
23 

9 
18 

9 
10 
15 
2 
20 

9 
18 
a 
27 
27 
37 
10 

6 
12 
13 
uy 
13 

7 
10 


1951 Jig 
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1933 XVI 
1937 II 24,0 
200 
220 
1924 II 40 265 
12 210 
8 235 
1950 I 235 25 225 <1,0 * bok 
1948 IV ESE 180 19 
ENE 125 12 24,5 
| 
1949 I ESE 230 15 240 <1.0 
ENE 200 10 245 
1948 X ESE 180 250 <1,0 
ENE 1,5 255 
1954 Edna ESE 260 250 (41.3 +2 
ENE 175 265 1.3 415 2 
1937 VIII SSE 130 275 <1,0 
ESE 230 250 
ENE 80 280 
1938 IV ESE 240 255 
ENE 185 260 
1944 I SSE 325 
ESE 310 275 
1951 Able ESE 255 275 <1,0 Toh 
280 
Teh 
ENE 65 310 
6.5 
ESE 195 "300 
| ESE 160 325 -15 7.0 \ 
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Octants NNE or NNW, a few of them were in ENE but far off from Charleston. 
Consequently none of them lowered sea level at Charleston appreciably. 
Tropical cyclone V 1927 caused a negative displacement of sea level of 
1.0 foot seven hours after a positive displacement of 1.9 feet. The negative 
displacement occurred while the cyclone center was in Octant NNW although 
it continued into NNE; however, the minimum distance to Charleston in NNE 
was 270 miles so the storm may have been out of effective range of Charles- 
ton by the time it reached NNE. Two other cyclones I 1934 and Able 1952 
passed through NNW also. Neither lowered sea level by a foot but V 1927 
was closer to Charleston while in Octant NNW than either of them. 


Tropical Cyclones in Octant WNW 


As was true of tropical cyclones in the previous category, tropical cyclones 
coming closest to Charleston in Octant WNW will ordinarily have Charleston 
on the right side of their track and therefore will cause positive displace- 
ments of sea level at Charleston because they have passed through Octant 
WSW and possibly SSW. However since they contribute little to raising sea 
level while in WNW, a negative time lag is to be expected. 

This is generally true of the cyclones listed in this category in the table. 
Four cyclones caused positive displacements of 2-1/4 to 4 feet with negative 
time lags of a few hours. Two others however, 1924 V and 1949 II, caused 
positive displacements of 2 feet with small positive time lags; these two 
were also the only cyclones in the group to cause subsequent negative sea 
level displacements. While the positive time lags may possibly be attributed 
to error introduced by inferring a constant speed between points of refer- 
ence, there is no apparent reason why only these two cyclones should cause 
negative displacements. Tropical cyclone V 1924 was very close to Charles- 
ton while in Octants NNW and NNE, but a number of the other storms were 
closer than Hurricane II 1949 while in these octants. Perhaps the unusual 
severity of Hurricane II 1949 accounts for its effect. 

The effect of Tropical Cyclone V 1924 was unusual too in that there was 
only a two hour difference in time between the maximum positive and nega- 
tive displacements of sea level. This is shown on the marigram as a re- 
versal of the tidal rise between 5 and 7 a.m. (75° W. time) on September 30. 
Immediately thereafter, the tide rose again at a somewhat accelerated rate. 
The eye of the cyclone almost passed over Charleston with strong winds first 
building up the level of the sea, the calm then releasing the accumulation and 
strong winds shortly thereafter driving the water from the coast. There was 
a change in sea level of 4.1 feet which took place in two hours that was direct- 
ly attributable to meteorological conditions. 

Hurricane XII 1933 caused a positive displacement of 1.0 foot with a time 
lag of -40 hours. The unusual time lag may be due in part to tide gage failure 
during the next eleven hours during which time the maximum recorded dis- 
placement may have been exceeded and in part to weakening of the hurricane 
as it approached Charleston; the storm dissipated shortly thereafter. 

One hurricane 300 miles from Charleston had little effect on sea level. 
The Charleston tide gage was inoperative during two other storms. 


: 
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Those cyclones coming closest to Charleston in Octant NNW have Charles- 
ton on the right side of their tracks. They may have large positive displace- 
ments of sea level with negative time lags based on their effect in previous 
octants. They may also have negative displacements with positive time lags 
based on subsequent effects in Octant NNE. 

Hurricane V 1928 caused a positive displacement of 3.8 feet followed by a 
negative displacement of 1.2 feet; the time lags corresponded to the above 
estimates. Hurricane I 1928 caused a displacement of +1.4 feet with a small 
negative time lag. 

Hurricane II 1929 caused a positive displacement followed by a negative 
displacement but the time lag for the rise in sea level is +4 hours. Here 
again the apparent discrepancy may be due to inference of a uniform speed 
between two points on the travel path twenty-four hours apart. The time lag 
for the negative displacement puts the cyclone center in Octant NNE at the 
time. 

Four other listed cyclones in this category had no appreciable effect on 
sea level at Charleston. Hurricane V 1941 is among these although it came 
within 20 miles of Charleston. Its failure to raise sea level is explained by 
its wide swing around Charleston in its clockwise motion but its failure to 
subsequently lower sea level while coming close to Charleston in Octants 
NNW, and NNE and ENE is not readily explained unless the hurricane lost its 
strength prior to approaching Charleston. 


Tropical Cyclones in Octant NNE 


Only one cyclone IV 1942, was closest to Charleston in Octant NNE. It 
was 240 miles away and had little effect on sea level at Charleston. 


Tropical Cyclones in Octant ENE 


Twelve tropical cyclones, all of hurricane intensity, were closest to 
Charleston in Octant ENE. Charleston was on the left side of the track of 
each. 

There was a positive displacement of between one and two feet by six of 
these hurricanes, each with a negative time lag indicating that they were set 
up while the cyclone was in the ESE direction. The negative displacement of 
1.5 feet for 1954 Hazel is interesting as it illustrates the behavior of a re- 
surgence as defined by Redfield and Miller. There was a positive displace- 
ment followed consecutively by a negative, a positive and then a negative dis- 
placement, each about six hours apart. These oscillations are believed due 
to the free motion of the water in restoring the disturbed levels, augmented 
in some cases by the changed direction of the wind. The period of oscillation 
following this storm, roughly that of the semidiurnal tidal cycle, appears to 
be the usual pattern at Charleston. If this is so the resurgence at Charleston 
is not apt to attain substantially greater heights than the initial rise as it will 
occur at about the same tidal phase. 
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Tropical Cyclones in Octant ESE 


August, 1957 


There are twenty-four tropical cyclones listed as coming closest to 
Charleston in Octant ESE. Of these only two have Charleston on the right 
side of their tracks. 

Although one of these raised sea level while the other lowered it, the data 
are not contradictory. Tropical cyclone III 1934 lowered sea level at Charles- 
ton by 1.2 feet. The negative time lag indicates it caused this while in Octant 
ENE. Hurricane V 1935 raised sea level 1.4 feet. The positive time lag in- 
dicates this occurred after the storm reached Octant SSE. 

Theoretically, tropical cyclones in Octant ESE will set up currents paral- 
led to the coast and have little effect on sea level at Charleston. This may be 
particularly true of cyclones with Charleston on the left side of their tracks. 
Some of the tropical cyclones listed in this classification caused positive dis- 
placements of sea level. In general these displacements are presumed to be 
due to their effect in octants previously traversed. The negative time lag or 
small positive time lag supports this suggestion. 

The displacements attributed to several tropical cyclones are probably due 
to extraneous conditions. Sea level rose 1.3 feet when tropical cyclone XVI 
1933 was near Charleston but this maximum displacement occurred about a 
half day after the cyclone had dissipated at its closest point to Charleston. 
There was a rise in sea level of 1.5 feet near the time Hurricane Jig 1951 
was first reported but equal displacements of sea level were reported on the 
previous two days. In both cases, the displacements appear to be coincidental 
to the presence of the tropical cyclones or at best they are only partially 
responsible. 

The two storms Hurricane II 1930 and Hurricane I 1945 illustrate the dif- 
ficulty in the classification of storm intensity. Both cyclones are listed as 
attaining hurricane intensity and traveled over almost identical paths while in 
the vicinity of Charleston. The first had little effect on sea level at Charles- 
ton and the second raised sea level 2.6 feet. Further investigation disclosed 
that the first had lost most of its strength before reaching the Atlantic and the 
second retained its hurricane strength, thus explaining the difference in effect. 


SUMMARY 


While it has been demonstrated that the bearing of a tropical cyclone and 
the direction in which the cyclone is moving with respect to Charleston are 
important factors in determining sea level displacement at Charleston, too 
many other variables are involved to permit accurate quantitative forecasts 
at present. It is possible to predict that certain storms will have little or no 
effect. It is also possible to predict that other storms will raise sea level by 
an estimated amount with the extreme coming within a few hours of a certain 
time, provided, of course that there is an accurate forecast of the cyclone 
track. 

As indicated previously, the intensity of a tropical cyclone at the time it 
can affect sea level at Charleston has not been adequately determined for the 
storms in this study. Similarly the dimensions of the effective diameters of 
most storms were not available although some data are published(2) and 
more information of this kind will result from the current National Hurricane 
Research Project. The speed of the cyclone centers in the various octants 
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are given in the table but no apparent relationship to sea level displacement 
has been noted. A tropical cyclone moving slowly should have a greater effect 
on sea level because winds will blow from the same direction for a longer 
period of time. On the other hand, a tropical cyclone moving with rapid 
velocity (forty to fifty miles per hour) will augment the velocity of winds on 
the right side of the track and thus may strongly affect sea level. Presuma- 
bly, heavy rainfall over the adjacent land area can raise the level of water in 
Charleston Harbor temporarily. 

All these factors really need to be considered to a certain extent in fore- 
casting sea level displacement. On the other hand, a study which classifies 
each cyclone according to each criteria may find that, with a limited sample, 
each storm is in a class by itself. In view of these variables, the fairly con- 
sistent results within each classification used in this study indicate the rela- 
tive importance of the criteria used. 
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ABSTRACT 


This paper describes the sand wave formations present in alluvial rivers 
and discusses their variations with respect to river stage. The related rise 
and fall of crossings is also discussed. The illustrations are prepared from 


fathometer surveys and show the changes in the sand waves due to changes in 
stage. 


INTRODUCTORY 


Close association with various aspects of planning, constructing and oper- 
ating works for the improvement of the Mississippi River, has gradually led 
to the conclusion that systematic changes take place in the river bed, and that 
these changes are a major factor in determining resistance to flow and in 
governing the general behavior of the river. Two somewhat over-lapping 
phenomena are involved; namely, (1) a rise and fall of the crossings between 
bends in sympathy with variations in discharge, and (2) a system of sand 
waves on the river bed, which varies as to general arrangement, and as to 
size and spacing of individual waves, in accordance with some controlling 
condition. 

As to the first phenomenon, the variation in height of crossings with re- 
spect to discharge, it will be noted that its existence is well known to practi- 
cal navigators, and engineers closely identified with the practical aspects of 
river work. Although obviously operating on all alluvial rivers, this phe- 
nomenon does not seem to be generally recognized among the hydraulicians, 
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and its effects have received little, if any, attention in evolving the theoreti- 
cal background of formal river hydraulics. 

As indicated, crossing bars, rise and fall in a manner broadly related to 
changes in discharge: On the lower Mississippi, the magnitude, of these 
changes ranges from 25%, to even as much as 50%, of the change in water 
surface level. This phenomenon is so sensitive that even a small rise is im- 
mediately reflected by a rise of the crossing bar. When a period of high flow 
is succeeded by a fairly rapid fall, these built-up crossings cause a river to 
be a series of pools with steep crossing slopes: Under the influence of the 
resulting high velocities, the crossings rapidly scour, with a resultant lower- 
ing of the pools to “normal” low water. Such scour may be disorderly, how- 
ever, and may fail to produce a good navigable channel; in fact, the more 
rapid the fall, the more likely it is that considerable dredging will be re- 
quired at the crossings. 

Closer study of this rise and fall of crossings, may throw new light on the 
movement of bed load material: It may develop that, in effect, the bends are 
“scoured out,” at high stages, and that much of the “scoured” material is 
simply moved to the first crossing. During lower stages, this material is 
probably gradually removed from the crossing and deposited in the bend be- 
low. In addition to this source of “back fill,” bends also receive back-fill 
ftom the extensive bank caving, which follows a falling river, after the sup- 
porting hydrostatic pressure is removed from the banks undermined at high 
stages. Consideration of these river bed “gyrations,” prompts a remark 
about hydrographic surveys; namely, that new bed configurations will be 
found just as often as such surveys are made, and the only really “complete” 
survey, would be a continuous motion picture of the river bed. 

This discussion of the rise and fall of crossings, in sympathy with changes 
in discharge, is presented for the purpose of relating it to the sand wave phe- 
nomenon, as one of the two major influences, affecting resistance to flow. It 
is hoped that the discussion will focus attention on a major hiatus in bed load 
Studies; that is, the special nature of bed load movement in bends, which is so 
unlike bed load movement over crossings (or in glass sided troughs in labora- 
tories), that it may be considered as virtually a completely different type of 
movement; that is, there may be an intermittent movement from bend to 
crossing (or vice versa depending upon stage) and a continuous “through 
freight” movement, going on in parallel. 

The conclusion cited, concerning sand waves and sand wave systems, was 
influenced by several empirical sources: namely, (1) long experience in 
hydrographic surveys, originally by piano-wire lead-line and circle -sounding 
equipment, and, since the 1930’s, by one of the early applications of sonic 
soundings; (2) routine sailing-line fathometer channel patrols for navigation, 
(3) evidence, from a steamboat pilot house, of uniformly spaced “reefs” on 
the convex side of bends, having (at this location) gradual upstream slopes, 
over which the water flows at high velocity, and short abrupt downstream 
slopes, forming comparatively deep and still pools, which sharply check the 
flow, and (4) experience in paddling a canoe, rowing a skiff, and swimming 
among such “reefs,” when all these means of locomotion gave personal con- 
firmation of the flow conditions described. 

The final conclusion is that these “reefs” extend into the channel, as deep- 
ly submerged sand waves, where they continue to produce the same effects as 
in the shoal water near the sand-bar shore; that is, they generate turbulence. 
This paper deals principally with these “reefs” as constituting a “primary” 
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wave system; that is, the “secondary,” and “tertiary” wave systems (both in 
evidence), are not given much consideration. The term “wave” will be re- 
stricted to these primary waves, and the secondary and tertiary waves will 
be referred to as “ripples.” 


Preliminary Field Investigations 


Continued speculation on the sand waves, led to the decision to make a 
preliminary field investigation: If the results warranted it, more extensive 
investigations could then be made, and any such preliminary work would 
determine the nature and direction of the more extended study. 

It should be noted here, that such investigations would be entirely imprac- 
tical except by the use of the fathometer to obtain a continuous profile along 
the path of travel; that is, no amount of conventional hydrographic surveying 
could yield the continuous detailed data necessary to study these sand waves. 
For this reason, previous “sand wave” studies, on the Mississippi, by 
McMath (around 1880) and Ockerson (about the turn of the century) are of 
little interest. 

A preliminary field investigation was arranged late in 1952: Advantage 
was taken of a regular Navigation Channel Patrol trip of the stern-wheel 
Steamer General John Newton to obtain a continuous fathometer profile of the 
“sailing-line” from New Orleans to Old River, a distance of some 200 miles. 
The horizontal scale ranged from 30 to 40 inches to the mile, and the vertical 
scale was one inch to 10 feet. The observations, on this trip, confirmed the 
general occurrence of sand waves over the crossings and around the bends. 


Hypothesis 


The following 4-point hypothesis was formulated on the basis of accumu- 
lated empirical experience and the results of the first fathometer “sailing- 
line” profile. 

Point (1) 

Sand waves occur generally on the bed of a river in a systematic manner, 
and they vary in a systematic manner with changes in discharge. 
Point (2) 

The sand waves tend to become larger with increases in discharge and 
smaller with decreases in discharge. 
Point (3) 

There is a lag in the rate of change of sand waves, and, on a continuously 
rising river, the waves are somewhat “smaller than they should be,” and, on 
a falling river, they are somewhat “larger than they should be.” 
Point (4) 

These waves constitute a major element of resistance to flow. 


Corollaries 


These corollaries proceed from the hypothesis: 


Corollary (1) 


Roughness, or resistance to flow, is not an accident or an attribute of the 
soil formations through which a river flows, but is rather a property, or 
condition, of the river bed, controlled by the river itself. 
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Corollary (2) 

The fact that the ‘rating curve,” at any given station, is a “loop” (Fig. 1) 
rather than a “curve” may be partially explained on the basis of a “lag” in 
sand wave adjustment; that is, the smaller waves present, on a rising phase, 
offer less resistance to flow, and the larger waves present, on a falling phase, 
offer more resistance to flow. 

The over-lapping phenomenon, of rise and fall of river crossings, also 
operates to produce the conditions recorded by the “loop” rating curve. 


Corollary (3) 

Slope itself (in a local sense) tends to become a “result” rather than a 
“cause.” 

Underlying all this discussion, is the assumption that, with discharge and 
other governing conditions constant, a stream will produce a system of sand 
waves, which remains constant, as to size and shape of individual waves, and 
as to the wave “system.” The implication is that, during a cycle of rise and 
fall, bed adjustment cannot take place, with sufficient rapidity, to be, at all 
times, in accord with the instantaneous demands of the varying discharge. 


Extended Field Observations 


Since the initial field observation in 1952, observations, on an extended 
basis, have fully confirmed that sand waves exist in a systematic manner, and 
vary in accordance with some orderly influence. Although the “systematic” 
occurrence is clear from visual inspection of the profiles, the influences 
governing the “system” are less clear. There are, however, definite indica- 
tions that “individual waves” and “wave systems” vary in sympathy with 
changes in discharge, and there is some likelihood that a certain type of wave 
and wave system will be found to be associated with the horizontal alignment 
at any location. 

As to the effect of increased discharge in increasing wave size, there was 
early evidence of an “anomaly;” that is, it was noted on two occasions, when 
the river rose to nearly bankful, that the “tops” of some crossings became 
perfectly smooth: The word “anomaly” is employed because waves were 
still present, on the upstream and downstream slopes of the crossings and 
around the bends, and the waves, where present, were larger than at low 
water. This behavior, of the tops of crossings, will be an interesting facet 
of study to pursue further, when high over-bank river stages again occur. 

As a first result of these field observations, the following tentative deduc- 
tions are advanced: The direct effect on resistance to flow, of the size and 
shape of the particles on the river bed (i.e., those on the immediate surface 
at a given instant), may be small in comparison with the effects of the suc- 
cessive orders of waves and ripples. Since the “primary” sand waves are so 
large, as to virtually constitute major river-bed “topography,” it appears 
possible that their frictional effect is great as compared with, either the ef- 
fect contributed by the lesser orders of waves (i.e., “ripples”), or the 
theoretical “drag” effect. Finally, it may be that the sizes and shapes, of the 
grains, are ultimately significant, not because of their direct effect upon fric- 
tion, but because of their indirect effect, in governing the formation, size and 
systematic variation of the primary sand waves and “ripples” on the river 
bed. It also appears likely that the size range of particles, present on the 
river-bed surface, at any given time, is subject to control by the river, itself, 
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as seems to be the trend of opinion in recent papers, based on laboratory and 
field measurements, and analytical studies. 

The observations, on which this paper is based, were entirely fathometer 
profiles; that is, no studies were made of sand sizes, the manner of adjust- 
ment of the waves, their speed of movement, or other associated phenomena. 


Views on Sources of Friction in Alluvial Rivers 


Various hydraulic texts and published papers, treating flow of rivers, bed 
load movement, and related questions, have been examined in order to 
develop the thinking with respect to systematic variations in the river bed 
(i.e., sand waves, and the rise and fall of crossings) as major phenomena 
governing resistance to flow and the general behavior of alluvial rivers. In 
general, the broad subject, of such systematic changes, does not seem to en- 
ter into formal river hydraulics: Rather, there is extended consideration of 
friction as such, and a change in attitude, over the years, as to the causative 
factors. 

The generally accepted view, of earlier writers, appears to have been that 
“roughness” is a property, inherent in the channel, and ultimately related to 
the nature of the soils, or the geological formations, through which the river 
carves its channel. The impression seems to have been that, whatever it was 
that controlled friction, it was constant for a given reach, or even an entire 
river, once its value had been determined. 

Another group of writers, relates the friction to the “sand paper” or “drag” 
effect of the particles, present on the surface of the bed. Some go further and 
indicate a belief that the river exercises a selective influence over the sizes 
of grains, actually present under any particular condition. Still other writers 
take a position, somewhat related to that presented in this paper; that is, sand 
waves and ripples, are a factor in determining resistance to flow. At least 
one has indicated the possibility that the waves get larger with increased dis- 
charge. Another thought the river bed would be smooth at the lowest stage. 

Still other writers make river bed behavior very complex: They put un- 
limited trust in involved mathematics which finally rests on field observa- 
tions, or laboratory measurements, that may contain an element of dubious- 
ness. By insisting on a detailed particularized approach, they implicitly 
reject the possibility of a generalized approach, based on laws or principles, 
perhaps susceptible of initial identification and subsequent measurement on 
an empirical basis. It may be, figuaratively speaking, that microscopic work 
on sand grains, and small areas of stream beds, has prevented recognition 
of (1) the need to view the river bed continuously in profile; (2) the possible 
significance of the systematic river bed changes discussed in the opening 
paragraph of this paper, and (3) the possibility that bed-load movement in the 
bends, is so different from bed-load movement in the crossings, that they 
constitute virtually separate and distinct phenomena. 

Summing up, it appears that the proposed hypothesis conflicts somewhat 
with both the older and the more recent prevailing views concerning causes 
of flow resistance, by implying that the major source of resistance arises 
from systematic river-bed changes, which the river makes, of its own voli- 
tion, and modifies to meet changing discharge requirements. The overall 
conception is a river which lets itself down its valley under control, much as 
a person controls his descent, down a flight of stairs, with some additional 
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control obtained by sliding his hand down the bannister. Final support, or 
partial refutation, of the four-point hypothesis, awaits additional fathometer 
observations at overbank stages, and comparisons of the resulting profiles 
with similar profiles, taken at lower stages. Because of the prolonged low 
stages, due to the drouth, no overbank stages have occurred in the five years, 
since interest in this subject came to a head. 


Comparisons of Fathometer Sailing-Line Profiles 


Successive sailing-line profiles were run, over several years, under vari- 
ous conditions of discharge. These profiles have been studied to detect sys- 
tematic variations, and to select half-mile “study-reaches” for comparing 
such variations. The results of these studies have been incorporated in the 
figures, now to be discussed: 

Each figure has a vicinity map to show the river configuration above and 
below the selected “study-reach.” There is also a condensed “sailing-line” 
profile, extending above and below the study-reach. The location of the study- 
reach itself is shown on the map and on the profile. Finally, the half-mile 
study reach itself is shown to full scale, for each of the conditions, to be 
compared. 

In viewing the full scale study-reach profiles, it should be noted that the 
sand waves are somewhat distorted: This is due to the fact that the “vertical” 
axis of the fathometer chart is a curve; hence, it is “vertical” only at the cen- 
ter of the chart, is inclined up and to the left on the lower half, and is inclined 
up and to the right in the upper half of the chart. The vertical scale is as in- 
dicated by the depth numerals, and the horizontal scale is as indicated by the 
marks, which set off 20ths of a mile, or roughly 250 feet. 

Figs. 2, 3, and 4 

A fairly typical crossing: In fig. 2, at the higher stage, the upstream slope 
has smoothed off, but there is a well developed wave system on the rest of the 
crossing, and the waves are larger at the high stage than at the low stage. 
Note that the amplitude of the waves, on figs. 3 and 4, is nearly 30 feet and 
that the normal primary waves have formed on these “Super” waves. It is 
suspected that the “Super” waves, just forming here, will take over complete- 
ly at higher than bankful stages. In fig. 4, the waves, on the downstream 
crossing slope, are of about the same amplitude at both stages. This devia- 


tion from the hypothesis, cited earlier, does not affect its general validity, as 
will appear in due course. 


Figs. 5 and 6 
Both of these figures deal with the lower end of a crossing at the start ofa 
bend. In fig. 5, the waves on the higher stage are larger than the waves on 


the lower stage. In fig. 6, the waves are much larger at the higher stage than 
at the lower stage. 


Figs. 7, 8, 9, and 10 

To the right of fig. 7, the Sand Wave System for the higher stage, has built 
up to an amplitude of about 10 feet. These waves are not present at the lower 
stage. Note that the crossing has smoothed off at the higher stage. In fig. 8, 
a “Super” Sand Wave System has attained an amplitude of over 20 feet, at the 
higher stage. Note that the normal primary waves are superimposed on the 
“Super” system. At the lower stage there is no significant wave system. In 
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fig. 9, the explanation for the system of large waves, at the lower stage, is 
that this is the same system of large waves, as shown in fig. 8, moved down- 
stream. As the stage drops, it is possible that the waves migrate down- 
stream, decreasing in size as they move. In fig. 10, the waves in these two 
systems are not too different; possibly because of the equalizing effect of the 
wave migration, as indicated in figs. 8 and 9. Note that the downstream 
crossing slope has been somewhat smoothed off, at the higher stage. 


Fig. 11 

This figure shows “Grandview Reach,” a long “reach” that has been his- 
torically stable. At the far left, for the highest stage, there is a well 
developed wave system with an amplitude of about 10 feet; at the middle 
stage, the wave system has migrated to the center of the study reach and the 
amplitude has decreased to about 4 feet; and, at the lowest stage, the major 


wave system has migrated further downstream and the amplitude has de- 
creased to 2 or 3 feet. 


Fig. 12 

This figure shows a rapidly caving bend above Baton Rouge: This bend is 
composed almost entirely of sand to the top of the bank. The unusual shape 
of the waves at the high stage is due to the reduced horizontal scale of this 
figure. The steep slope from the lower end of the crossing to the head of the 
bend has smoothed off at the high stage. This bend has caved at the average 
rate of 5,000,000 cubic yards per year for the last 16 years and has moved 
back 70 feet per year at the point of greatest caving. 


Figs. 13 and 14 
Note the unusual alignment of the river in this vicinity. The waves are 


larger on the higher stage. In fig. 14, there are indications of “Super” waves, 
on which the normal primary waves are superimposed. 


Figs. 15 and 16 

These figures show a typical crossing: Fig. 15 shows the upstream cross- 
ing slope, and fig. 16 the downstream crossing slope. In fig. 15, there isa 
well developed “Super” wave system, on which are superimposed the normal 
primary waves. However, at the lower stage, only the normal primary waves 
are present. Fig. 16, shows rather symmetrical and consistent wave systems, 
at both stages. There are also indications that there has been an “equaliza- 
tion” of the wave systems, at higher and lower stages, due to the migration of 
the larger waves from the upstream crossing slope. 


CONCLUSIONS 


These studies afford substantial support for three of the component points 
of the hypothesis; namely, (1) that sand waves exist generally on the river bed 
in a systematic manner; (2) that sand waves tend to vary in amplitude and 
spacing, in accordance with variations in discharge; and (3) that a major 
source of resistance to flow, in an alluvial river, is provided by a primary 
system of sand waves, which, in turn, is created by the river itself and oper- 
ates as a means of controlling its orderly progression down its valley. The 
fourth “point,” that the waves constitute a major component of resistance to 
flow, has not been confirmed, although it still appears reasonable. 

The fathometer profile observations, on which this paper is based, covered 
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a limited range of river stages, and the value of the conclusions, derived 
therefrom, is correspondingly limited. The paper will serve its purpose if it 
demonstrates the existence of a practical field of investigation that may lead 
to useful knowledge concerning; (1) the adjustment of sand wave formations, 
(2) the rise and fall of crossings, (3) the special character of bed load move- 
ment in bends, and (4) the joint bearing which these phenomena have upon re- 
sistance to flow. 

Federal Agencies, engaging in planning, constructing and operating water 
resources projects, can often gather data on river bed conditions at little ex- 
pense over that of the fathometer paper involved; as, for instance, in the case 
of channel patrols for navigation. In general, observations can be made at 
full cruising speeds on trips for other purposes, or on the idle return trips 
from other missions. Either the “sailing-line,” the “thalweg,” or the river 
“center-line,” can be followed, as may be convenient, and it is not necessary 
to cover an “area” in order for the data to be useful. Working up the data, 
can conceivably be relatively inexpensive. A point to note is that studying 
“lines” or “profiles,” has much to be said for it, as against studying “cross- 
sections”: In the first place, the key to studying these river bed phenomena 
is “continuity,” and this is something “cross-sections” cannot supply, despite 
the great amount of field and office work involved; in the second place, a 
“line” can be taken as representative of a “volume” within certain limits, if 
proper lateral-control of the line is provided, and known points are selected 
for linear-control “fixes.” 

After a certain amount of “qualitative” study of “lines,” to determine the 
manner and magnitude of river bed adjustments, ways and means will prob- 
ably be evolved to use the same data to make “quantitative” determination of 
the bed-load material movement involved in the sand waves and in the rise 
and fall of crossings; that is, the material that probably has the most pro- 
found effect on the hydraulic properties of streams. 

In conclusion, it appears likely that organized study of these phenomena of 
river bed behavior, might have these effects: (1) throw new light on some 
theoretical aspects of river hydraulics, (2) indicate a better approach to plan- 
ning and constructing bridge piers, pipe-line crossings, etc., and (3) point the 
way to a better understanding of the hydraulic results of comprehensive river 
stabilization and lay the basis for extending the claimed economic benefits, 
for such work, far beyond the mere protection of lands and improvements 
from the effects of bank caving and bar building. 
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HYDRAULICS DIVISION 


SYNTHETIC STORM PATTERN FOR DRAINAGE DESIGN 


Clint J. Keifer! and Henry Hsien Chu,” Associate Members ASCE 
(Proc. Paper 1332) 


SYNOPSIS 


The City of Chicago has been engaged in an extensive auxiliary outlet 
sewer construction program since the end of World War II. Bond issues 
totaling more than $88,000,000.00 have been passed to finance this program. 
The engineers of the Sewer Planning Division felt it worthwhile to make a 
comprehensive study of the principles involved in the rainfall-runoff relation- 
ship in order to economically utilize this expenditure. The development of a 
synthetic storm pattern, for use in the hydrograph method of sewer design, is 
one phase of this research, and is the theme of this paper. 


INTRODUCTION 


It is the purpose of this paper to present a method of determining a storm 
pattern which may be used in the hydrograph method of sewer design.(1,2,3) 
This storm pattern will be developed from a rate-duration-frequency curve 
selected for use in design. It will also be shown that although this storm pat- 
tern will include the average intensities of the rate-duration curve for all 
durations, it will produce no greater runoff peaks than that of the rainfalls 
of the separate durations. 

The chronological location of the peak period of rainfalls with reference to 
the total storms period and the amount of antecedent precipitation immediate- 
ly preceding the maximum period of any durations will be derived from the 
Statistical average of rainfall records. The importance of this antecedent 
rainfall in affecting the peak of the runoff hydrographs will be illustrated. 


Note: Discussion open until January 1, 1958. Paper 1332 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 4, August, 1957. 


1. Senior Sewer Designing Engr., City Bureau of Eng., Chicago, Ill. 
2. Civ. Engr., City Dept. of Public Works, Chicago, Ill. 
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The rainfall and resulting runoff hydrograph should be analyzed in order 
to determine what other characteristics, besides the average rate of water 
falling within a given time period, would tend to affect the peak rate of run- 
off. If in one rainfall, some precipitation preceeded the maximum period, 
filling many of the surface depressions, wetting all surfaces and satisfying 
the greater infiltration capacity of the soil, the peak runoff rate would be 
higher than from another rainfall having the same intensity during the maxi- 
mum period but lacking the antecedent rainfall. In the development of a syn- 
thetic storm pattern for use in the hydrograph method, the right amount of 
antecedent rainfall must be replaced in front of the statistical volume of rain- 
fall occurring within the maximum period. 

Another important characteristic that affects the peak rate of runoff with a 
given volume of water falling in a specific period, is that of the variation of 
rainfall intensities within that time period. Consider a long duration rainfall 
in which most of the precipitation falls in the early part of the period. In this 
case, the greater infiltration capacity at the beginning and the surface depres- 
sion storage would absorb much of the peak rainfall, thus making the runoff 
rate approach a lower value of equalibrium. On the other hand, if the bulk of 
precipitation occurred in the latter portion of the period, most of the previ- 
ously mentioned losses would be already satisfied before the time of peak 
rainfall intensity, and a higher peak rate of runoff would result. 

To sum up briefly, the three most important characteristics affecting the 
peak runoff rate for a specific period or duration are as follows: 


1. Volume of water falling within the maximum period, 
2. Amount of antecedent rainfall, 
3. Location of the peak rainfall intensity. 


A synthetic rainfall encompassing the statistical average of these three char 
acteristics should then be adequate for use in the hydrograph method as the 
design storm pattern for any given duration, see Figure (1). 


Evaluation of the First Characteristic for Single Storm Pattern 


The volume of water falling within the maximum period can be taken from 
the rate-duration curve of a given frequency selected for design purposes. 
The equation for this rate-duration curve may take the following form: 


(1) 


where igy, is the average intensity in inches per hour, tg is the duration of the 
maximum period in minutes, and a, b and c are constants. The volume of 

water falling, or mass rainfall for any duration is graphically represented by 
the area of a rectangle of which the upper right-hand corner touches the rate- 


duration curve, see Figure (2). This rainfall mass, in inches, can then be ex- 
pressed by 
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First, a completely advanced type rainfall pattern can be developed which has 
the peak occurring at the beginning, and with no antecedent rainfall. A curve 
representing the hyetograph of this type of rainfall must have the same aver- 
age intensity for any maximum period as that given by the rate-duration 
curve, as shown in Figure (2). The area under the curve of the hyetograph 
can be expressed mathematically by 


4 


where i is the ordinate of the hyetograph in inches per hour. 
Differentiating with respect to tg 


But from Equation (2) 


Differentiating 


dt, 60 (4) 


Combining Equations (3) and (4) 


a (1-6) t c] 
(5) 
Equation (5) then represents a completely advanced type storm having for all 
durations, the same average intensities as those given by the rate-duration 
curve. For Chicago area, the accepted rate-duration curve developed by 
Eltinge and Towne of the Sanitary District of Chicago, (4) for a five-year fre- 
quency, is 


90 


(6) 


Therefore, the constants in Equation (1) are a = 90, b = 0.9, and = 11. 


dP = (3) 
dt, 60 
60 
. 
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Substituting these values into Equation (5), one obtains 


90 (0.1 (7) 


By assuming different values for tg, Equation (7) represents a storm starting 
with the rate of 8.18 inches per hour at the first instant and decreasing 
steeply to 4.42 inches per hour at the end of the 5th minute, to 1.13 inches 
per hour at the end of the first half hour and so on. 

To make Equation (5) applicable to any intermediate type storm pattern 
the following modifications shall be made: Within the maximum period of 
any rainfall, the duration tg can be split up into that part occurring before 
the most intense moment and that part after the most intense moment. In- 
troduce the symbol ‘r” to represent that portion of any duration occurring 
before the most intense moment, expressed as a ratio to the entire duration, 
refer to Figure (3), and set 


rt and = (l-r)t (8) 


where tp is the time before the peak in minutes, and measured from the peak 
to the left and t, is the time after the peak in minutes measured to the right 
of the peak. The ratio ‘r” is then a measure of the advanceness of the storm 
pattern. If r = 0, the storm is of the completely advanced type as discussed 
in the last paragraph. If r = 1, the storm is of the completely delayed type, 
that is, it has its peak at the end of each and every duration and has con- 
siderable antecedent rain falling before each and every maximum period. 

For intermediate type storm patterns, the “r” value will range between zero 


and one. 
Solving the first part of Equation (8) for tg and substituting in Equation (5) 
-b)(% 
(Before The Peak) a((! Me) c] (9) 
+ 
Solving the second part of Equation (8) for tq and substituting in Equation 
(5) 
b 
+c] 
(After The Peak) ‘= a (10) 
¢] 
A synthetic hyetograph plotted from Equations (9) and (10) will have, for 
all durations taken during the most intense period, the same average intensity 
as the rate-duration curve from which the constants a, b and c are derived. 
Using the above mentioned constants for the Chicago area, 
t 0.9 
(Before The Peak) ) (11) 
r 
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t, \o9 
(After The Peak) 90 [0.1 ul] (12) 


t, \o9 2 
Determination of “r” Conforming with the Second and Third Characteristics 


The first characteristic has been evaluated by the development of Equa- 
tions (11) and (12) from the accepted rate-duration curve. It remains to de- 
termine, from records of excessive rainfalls, the best value to assign the 
constant ‘r” so that the other two characteristics of the synthetic storm pat- 
tern will conform approximately to the statistical data. 

Rainfall records were obtained from the Metropolitan Sanitary District of 
Greater Chicago, where the records of heavy storms from seventeen automa- 
tic recording rainfall gauges located in the Chicago area are compiled. 

Four widely separated stations were selected so as to retain the individuality 
of each station-rainfall. The four stations were at Loyola University (Station 
10), Springfield Avenue Pumping Station (Station 14), The West Side Treatment 
Works (Station 7) and Roseland Pumping Station (Station 15). Only those 
rainfalls which at least for some duration are considered excessive, were re- 
corded. Eighty-three station-rainfalls are given, with data taken from the 
records for durations of 15, 30, 60, and 120 minutes. 

In Table 1, the following information is given: For example, at Station 14 
on May 11-12, 1935, the maximum rainfall occurring in a fifteen minute 
period was 0.39 of an inch, while the antecedent rainfall occurring before that 
period was 1.40 inches. Since the maximum fifteen minute period contains 
three five-minute periods, the most intense was the second. Likewise, in the 
same storm at Station 14, the maximum mass occurring in the 120 minute 
duration was 1.82 inches, only 0.05 inches fell antecedent, and out of the 
120 minute duration the twentieth five-minute period was the most intense, 
therefore the storm was of the delayed type pattern. The blanked spaces indi- 
cate those durations within which the maximum mass rainfall was not exces- 
sive. These figures are purposely deleted in order that the statistical aver- 
age obtained will reflect the characteristics of only the excessive storms 
concerned. The mean values of antecedent rainfall and the location of the 
peak for each of the given durations were obtained. The mean values of mass 
antecedent rainfall are plotted in Figure (4), and the mean values for the loca- 
tions of peak are plotted in Figure (5). 

An analytical expression for the antecedent rainfall of the synthetic storm 
pattern can be derived as follows: For a given ‘r” value, the mass rainfall 


within the duration tq occurring before the peak can be obtained by multiply- 
ing Equation (2) by “r,” or 


la 


where Py is the mass rainfall before the peak in inches. 
For the Chicago five-year frequency rainfall 
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Then the total duration of the synthetic storm pattern is set equal to the long- 
est time of concentration of any sewer system to be designed. For Chicago 
the total duration of the synthetic storm pattern was set at 180 minutes. Also, 
since this duration was quite long for storm sewer design, it was assumed 
that any small amount of rainfall preceding this period would produce little 
effect on the rate of runoff, therefore no antecedent rain was assumed before 
this 180-minutes-long synthetic storm pattern; then the amount of antecedent 
rainfall preceding the duration tg, can be expressed as 


A=(R] -[R] = = (13) 


where A is the mass antecedent rainfall in inches. Solving Equation (13) for 
“r” and substituting the statistical mean values of antecedent rainfall for A 


and the corresponding duration tg, the value of “r” for each duration can be 
computed, see Table I. 


| 


50.8 | 1.769 


+ 


2.110 -180 


“ 
Table II = Determination of the Value of fr from Mass 
Antecedent Rainfall. 


Since the “r” is to be constant through all durations, and the antecedent rain- 
fall for the shorter durations is much greater and of greater importance in 
affecting runoff than that for the longer durations, it seems that the “r” com- 
puted for different durations should be weighted in proportion to the mass it 
reflects in obtaining the average ‘r” value. The weighted average value of 
“r” thus obtained is 0.386. 

The location of peak within each duration provides another direct approach 
in evaluating the “r” value. In the last line of Table I, “Mean Location of Peak 
within Rainfall Duration” or mean “r” for each duration, is calculated on the 
basis that the momentary peak of every recorded rainfall happened exactly at 
the middle of the recorded “Peak Five-Minute.” Then these “r” values of 
each duration were weighted proportionally to the time of duration in obtaining 
the average “r” for all durations. This was done because the shorter dura- 
tions contain so few five-minute periods that the mean location of the peak is 
not so well defined as in the cases of longer durations. The value of ‘r” thus 
derived by considering the location of the peak is 


| 
9 9 | A | | 
t | 4 A ¢ | Ar 
d d d ty | r | } 
15 | 11.5 | 22.5 | 1.000 1.290 | .416 | .322 | .1340 
30 21.4 | 32.4 | 1.388 | .902 | .324 | | 
120 | 74.4 | 65.4 | § .118 | .655 | .0771 
A=1.064 | )=.4115 
Weighted Value of “r”=.386 
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(.453) + 30(4/5) + 60(.356) + 120(. 366) 


75+ 304 G0 F120 


An even fraction 3/8 or 0.375 which is close enough to the above findings 
was assigned to “r,” thus establishing a synthetic storm pattern with a three- 
eighths advanced type. The mass antecedent rainfall for “r” = 0.375 is plot- 
ted in Figure (4) by the use of Equation (13). Also in Figure (5) a horizontal 
line with “r” = 0.375 is drawn. Comparison of these curves with the statisti- 
cal mean values indicate fairly close agreement in antecedent rainfall with 
the shorter durations and in the location of peak with the longer durations. 
Substituting the value of “r” = 3/8 into Equations (11) and (12) they become 


(8 0.9 
(After the peak) + 


The synthetic storm pattern plotted according to the above equations satis- 
fies best the three characteristics outlined in the beginning of this treatise 
and therefore can be used as the design storm pattern for the hydrograph 
method of storm sewer design in the Chicago area. Since the total duration 
of the synthetic storm pattern is limited to 180 minutes, the time before the 
peak is limited to 3/8 of 180 minutes or 67 1/2 minutes and the time after the 
peak to 5/8 of 180 minutes or 112 1/2 minutes. A hyetograph of the synthetic 
storm pattern is shown in Figure (6). 

Also shown for comparison, in Figure (6), is a hyetograph for a one-year 
frequency synthetic rainfall, based upon the Eltinge-Towne formula for a 
“One-year frequency” rate-duration curve, namely 


and the value of “r” is assumed to remain at 0.375. 


CONCLUSION 


The “Synthetic Storm Pattern,” having been derived from the accepted 
rate-duration curve, and fitted to the statistical values of antecedent rainfall 
and location of the peak, may now be used in the “Hydrograph Method” of 
sewer design. The advanced type pattern with “r” = 3/8, found to be accept- 
able in the Chicago area, might not be suitable in other localities. The au- 
thors believe that setting the “r” value is of considerable importance since 

it determines the amount of antecedent rainfall for all durations. In Figure 
(7), it is shown for comparison, the affect of the change in the “r” value on 
the runoff hydrograph of a 70 foot length of turf strip after subtracting losses 
due to infiltration, depression storage and surface detention. It can be seen 
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that for the completely advanced storm pattern, the peak runoff rate is only a 
third of that for the 3/8 advanced storm pattern. Figure (8) shows a compari- 
son of the runoff hydrographs from a typical residential area in Chicago, 
with 36 percent imperviousness, resulting from synthetic storm patterns hav- 
ing “r” = 0.375 and *r* = 0. 

The Synthetic Storm Pattern having the same average intensities as that 
given by the rate-duration curve for all durations, would not be likely to occur 
at the same frequency as the rate-duration curve from which it was derived. 
It would undoubtedly have a greater return period, that is, less frequent occur- 
rence. The use of the Synthetic Storm Pattern, however, should not produce 
greater runoff peaks then the use of a separate uniform intensity rainfall for 
every location in the sewer system, if the uniform intensity rainfalls have the 
correct antecedent precipitation. In Figure (9), the synthetic storm pattern 
has been used as a supply hydrograph at intervals along a sewer, and routed 
by the time-offset procedure to produce hydrographs at various times of 
travel. This use of the synthetic storm pattern, and subsequent use of uni- 
form rainfalls as the supply hydrographs, is only illustrative and would occur 
only if the area was entirely impervious and had no losses or reductions be- 
fore entering the sewer. Figures (10) and (11) show various duration uniform 
intensity rainfalls, with antecedent precipitation equal to that of the 3/8 ad- 
vanced synthetic storm pattern, being routed by time-offset to produce hydro- 
graphs at various times of travel. Figure (12) shows the comparison of the 
peaks of the hydrograph derived from the synthetic storm pattern and those 
derived from the uniform intensity rainfalls. It can be seen that the attenua- 
tion of the hydrograph peaks of the synthetic storm pattern envelope those of 
the uniform rainfalls of separate durations, and would produce no greater 
quantities for use in design. In a more rigorous analysis of comparing the 
design rates derived from the synthetic storm pattern to those derived from 
rainfalls of separate durations, it was found that the peak rate of runoff is 
directly proportional to the net volume of the supply curve in the maximum 
period equal to the routing time. Then, since the synthetic storm pattern 
has the same volume of rainfall up to the end of any duration as the uniform 
intensity rainfall of the same duration with antecedent precipitation, they 
would both produce the same runoff peaks. If infiltration, depression storage 
and surface detention were subtracted chronologically from these two rain- 
falls, the same net volume of supply would result. Thus after routing, they 
would again produce equal runoff rates. 

From the foregone, it is concluded that the synthetic storm pattern can be 


used as the “Design Storm Pattern” in the hydrograph method of sewer design 
for all locations in a sewer system. 
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ABSTRACT 


Generalized experimental results are presented to show the variation in 
both the discharge coefficient and the nappe profile for two-dimensional flow 
over a vertical sharp-crested weir at the end of a horizontal channel as the 

ratio of head to depth of flow changes continuously from Zero to unity. 


SYNOPSIS 


Whereas most existing weir formulas lose significance as the ratio of head 
to height of weir becomes great, it is known that the rate of flow over a sill at 
the end of a channel can be estimated from the basic weir equation by assum- 
ing the velocity of approach to have its critical magnitude. This gives rise to 
a discharge function for sills which meets that for weirs in a sharp peak when 
the head is about ten times the weir height. There are related indications 
that the nappe attains a maximum elevation and thickness at the same head- 
height ratio. Laboratory tests to determine quantitatively the variation in 
discharge coefficient and nappe profile over the entire range are described, 


and the experimental results are presented in generalized form for design 
purposes. 


Historical Background 


Poleni,(1) in 1717, was the first to publish an equation for weir discharge 
obtained by integrating over the outlet section an expression for the velocity 
in terms of the square root of distance below the free surface. Although he 


Note: Discussion open until January 1, 1958. Paper 1345 is part of the copyrighted 
Journal of the Hydraulic Engineering Division of the American Society of Civil 
Engineers, Vol. 83, No. HY 4, August, 1957. 


1. Iowa Inst. of Hydr. Research, State Univ. of Iowa. 


2. Prof., Fluid Mechanics & Div., Inst. of Hydr. Research, State Univ. of 
Iowa. 
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neither took account of the nappe contraction nor explicitly introduced the 
gravitational acceleration g, Poleni’s name continues to be associated with 
the elementary discharge equation 


in which q is discharge per unit width of flow section, C, the coefficient of 
contraction, and h the head on the crest. Contraction of the nappe was taken 
indirectly into consideration by Du Buat(2) sixty years later in setting as ar- 
bitrary limits of integration the crest level and a level the distance h/2 above 
the crest; it was also Du Buat who introduced present-day methods of estimat- 
ing the effect of submergence. Bidone,(3) better known for his study of the 
hydraulic jump, appears to have been the first to show an interest in express- 
ing algebraically the geometry of the weir nappe. 

Not until about the middle of the 19th century was attention focused upon 
the velocity-of-approach “correction” for weir discharge, which was to domi- 
nate certain branches of the literature for another century. In 1845 
Weisbach(4) published the familiar integral form 


eg 2g 


Since the solution of this equation has to proceed by successive approximation 
because of the implicit presence of the approach velocity vg on the left side 
of the equation, the impression has persisted that this velocity is independent- 
ly variable. Weisbach himself appreciated the fact that a unique solution 
existed for any ratio of head h to depth of approach h + w, as he proposed an 
empirical formula for its evaluation: 


2 
2 h (3) 


Bazin(5) in 1888 adopted the same form of function, with specific evaluation of 
the limiting coefficient of contraction (Cy)o 


4 
q= (0.405 95008) + 0.55 5 13/2 (4) 


the dimensional term 0.003/h being introduced to compensate for low-head ef- 
fects now attributed to surface tension and viscosity. In a series of investiga- 
tions extending over the next ten years, Bazin also went deeply into such 
varied aspects of weir flow as the distribution of velocity and pressure 
through the nappe, and the effects of incomplete ventilation, submergence, 
weir inclination, and geometry of the crest; his nappe-profile measurements 
are still used as the basis of spillway design. 

Bazin’s formulation of the discharge coefficient soon led others to attempt 
further improvements by slight changes in the type of algebraic expression 
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used. From 1911 to 1928 Rehbock(6) proposed a number of variants, the 
most pertinent of which was 


a = (0.605 + 0.03 2 + 2 (5) 


Comparison of Eq. (5) with Eq. (4) will show that the primary distinction lies 
in the use of the ratio h/w in place of h/(h+w) as the geometric parameter. 
The remaining term of Eq. (5)—which, like the corresponding one of Bazin, 
is not dimensionally homogeneous with the others—again reflects capillary 
and viscous effects at low heads. If the latter factors are ignored, and if the 
basic weir equation is written in terms of a discharge coefficient Cq, 


6 
a= (6) 


it will be seen that the essential aspects of the Bazin and Rehbock formulas 
may be reduced for comparison to 


= 0.608 + 0.384 (2) (7) 
h+w 


C. = 0.605 + 0.08 2 


a (8) 


which are plotted against h/(h+w) in Fig. 1. The two curves evidently differ 
by little more than 1% for heads as great as twice the height of weir. 

When Mises,(7) in 1917, determined the coefficient of contraction analy- 
tically for various two-dimensional forms of conduit outlet, he noted that the 
values for a symmetrical slot in a normal plate at the end of a conduit would— 
in combination with Eq. (2)—duplicate the trend of the Rehbock formula very 
closely over a considerable range if the ratio of slot dimension b to conduit 
dimension B were used to represent the ratio h/(h+w). His curves for Ce and 
Cg are also shown in Fig. 1. To be noted in particular is the limiting the- 
oretical value 7/(7+2) = 0.611 in comparison with the empirical values 0.608 
and 0.605 chosen by Bazin and Rehbock. One of Mises’ students, Lauck,(8) 
nearly a decade later approximated not only the discharge coefficient but also 
the entire nappe profile for the weir of infinite relative height, by means of 
the Cauchy integral theorem. 

Inspection of the Bazin and the Rehbock formulas will indicate that the lat- 
ter must approach the limit Cg =o as the relative height of the weir ap- 
proaches the limit zero (i.e., the free overfall), a trend that is necessarily 
shared by the orifice function of Mises, whereas the Bazin function has the 
finite limit Cg = 0.942. Now it is obvious that discharge over the free over- 
fall is finite rather than infinite, since the effective head is the critical depth. 
However, it was noted by Boss(9) that the critical stage of flow is already es- 
tablished upstream from terminal sills well before they attain zero height, 
which would yield a trend contrary to that of Eq. (6). That is, following the 
procedure of Rouse,(10) the discharge coefficient for such conditions can be 


and 
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determined by introducing into Eq. (6) the critical-depth relationship 


z 9 
= 


(10) 


which is plotted at the upper right of Fig. 1. Evidently, the discharge char- 
acteristics of terminal sills follow quite a different functional relationship 
from that of even moderately low weirs, with an intermediate maximum value 
of the discharge coefficient. Now the coefficient of contraction C, can be de- 
termined from simultaneous solution of Eqs. (10), (6), and (2) and plotted as 
at the lower right of Fig. 1. If, after Rouse, this coefficient is assumed to 
approximate the relative thickness of the nappe at its highest point, then this 
plot indicates that the flow profile will also attain a maximum stage at the 
borderline between weirs and sills. 

In order to investigate these aspects of the problem, the second author 
undertook in 1931-33 in the River Hydraulic Laboratory of the Department of 
Civil and Sanitary Engineering at M.I.T. a systematic series of profile (as 
well as velocity and pressure) measurements on the nappe of a weir that was 
successively reduced in height until the overfall limit was reached. The re- 
sults of the profile tests were reflected in a subsequent book,(11) but they 
were found to cover the transition range inadequately and hence were never 
published in their entirety. Similar measurements over a less extensive 
range of the parameter h/(h+w) were undertaken later by the Bureau of Re- 
clamation and published(12) for purposes of spillway design. Because it was 
believed that certain characteristics of the weir-sill range were nevertheless 
of technical interest, generalized plots of the M.I.T. data showing a parallel 
between the discharge-coefficient and nappe-profile functions were included 
in the Proceedings of the Fourth Hydraulics Conference(13) of the Iowa Insti- 
tute of Hydraulic Research. 

Continued efforts have since been made at the Iowa Institute to complete 
both the analytical and the experimental phases of the weir and sill problem. 
The relaxation method has already been applied to the improvement of 
Lauck’s approximation of the nappe profile for the infinitely high weir(14) and 
at the time of writing the study is being extended to other values of h/(h+w). 
However, since numerical integration is not as readily applicable to low weirs 
and sills as to high, the gaps left in the M.I.T. study were in 1956 filled in ex- 
perimentally by the first author(15) under the direction of the second, to the 
end of providing a basis for both design and further theoretical analysis. It 
is to the presentation of the composite results that this paper is devoted. 


Experimental Procedure 


The M.I.T. experiments were conducted in a 50-cm glass-walled flume, 
the rate of flow through which was measured with a ventilated sharp-crested 
weir 40 cm high in accordance with Eq. (5). The experimental weir was 


Thus, 
3/2 
= 1.06(1 +¥) | | 
ray. 
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initially of identical characteristics, but its effective height was changed by 
the addition of successive false floors on the upstream side. These were 
constructed of troweled cement mortar over sand, and each consisted of a 6- 
ft. level section preceded by a very gradual rise from the original floor of 

the flume. Three weir heights—40, 20, and 10 cm—were produced in this 
way. The 40-cm weir was then removed, and still another false floor was 
added at a height of 37.5 cm above the flume floor, with a brass-angle end to 
simulate the free overfall. A supplementary weir crest could be attached to 
the end with wingscrews to produce additional weir heights of 5, 2 1/2, and 
11/4 cm. Full ventilation was provided below the nappe in all instances. 

Profile measurements were made at 1-cm sections upstream and down- 
stream from the crest with point and hook gages reading by vernier to 0.1 
mm. For each of the seven weir heights the same nine rates of discharge 
were investigated; these corresponded to heads on the measuring weir rang- 
ing by 2-cm increments from 2 to 18 cm. Supplementary measurements of 
pressure distribution along the floor, up the weir plate, and through the nappe, 
as well as checks upon the constancy of the total head at various points, were 
also made, but neither the experimental techniques nor the results are per- 
tinent to the present discussion. 

The tests conducted at the Iowa Institute employed the same type of free- 
overfall structure as the one finally constructed at M.I.T., except that it was 
built of brass plate and angles and mounted in a 1-ft. flume. The level por- 
tion was 4 ft. in length, its height above the floor of the flume was 1 ft., and 
the approach curve had the form of a quarter ellipse with an axis ratio of 
3.5:1. The adjustable weir plate was machined to a sharp upstream edge, the 
top surface having a width of about 1/32 in. To minimize capillary distortion 
of the lower surface of the nappe at small weir heights, the crest surfaces 
were rubbed with paraffine. For the very smallest height (0.005 ft.) the crest 
was remachined to a knife edge. .The test procedure was essentially the same 
as that at M.I.T., although the range of head was limited to 0.1 ft.ch<0.5 ft. 
and the range of weir height to 0.005 ft<w<0.05 ft. Discharge measurement 
was effected by means of a 90° elbow meter in the 6-in. supply line that had 
been previously calibrated against a 12-in. sharp-crested weir mounted in 
the flume itself, the Rehbock formula being used as reference just as in the 
M.I.T. experiments. 


Discussion of Results 


As has long been realized, the head on a weir must be measured a con- 
siderable distance upstream from the crest, where the effect of the curvi- 
linearity of flow upon the surface elevation is still negligible. A distance 
equal to three or four times the head is generally considered to suffice. Un- 
fortunately, with relatively low weirs the effect of channel resistance becomes 
appreciable, so that the nappe profile is asymptotic to a sloping rather than a 
horizontal line. In other words, the head measurement then is subject to two 
types of error, the one increasing with proximity to the crest and the other 
with distance from it. 

In order to determine that section at which the net error would be a mini- 
mum, so that the heads for different weir heights could be evaluated con- 
sistently, a composite dimensionless plot was made of all profile measure- 
ments upstream from the crest for the most extreme case—the free overfall— 
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as shown in Fig. 2. Since the computed critical depth was used as the refer- 
ence length, and since the critical section has been shown to be the last at 
which hydrostatic distribution of pressure (i.e., absence of curvilinearity) 
can be assumed to prevail,(10) the abscissa of the desired section should be 
that having an ordinate value of unity, or x/h = 4. Inasmuch as this value for 
the free overfall agrees with that commonly used for weirs of great relative 
height, it was adopted for all subsequent evaluations in the intermediate range 
as well. 

From the experimental determinations of discharge and surface profile for 
the various heights of weir and sill, the coefficient of discharge was evaluated 
according to Eq. (5) and plotted as in Fig. 3. Herein the Rehbock parameter 
h/w is used as abscissa in the weir range, because of the linearity of the re- 
sulting function. Its reciprocal w/h then logically serves as the abscissa in 
the sill range, the relative scales of the two being so chosen as to yield ap- 
proximate agreement in the intermediate range. All of the Iowa data are 
shown, together with a comparable number from M.I.T. (i.e., all at the higher 
heads). Inspection of the diagram will indicate that Eq. (8) for weirs applies 
with good approximation over the range 0<h/w<6, and Eq. (9) for sills over 
the range 0<w/h<0.06. Although some discrepancy between the two sets in 
the zone of overlap is apparent, a smooth transition between the two functions 
can readily be drawn. 

In Fig. 4 are superposed a systematic series of nappe profiles extending 
through the weir range to the case of maximum coefficient, h/w = 10. Figure 
5 shows the corresponding family of profiles for sills. Like the coefficient 
plot of Fig. 3, the profiles (except for the lower surface of the weir nappe) 
are seen to be characterized by similar trends, the case h/w = 10 or w/h = 
0.1 evidently representing a maximum for both the discharge coefficient and 
the nappe deflection. 

As had previously been noted by the second author and used as the basis of 
the generalized profile diagrams already mentioned,(13) plots of z/h as func- 
tions of h/w and w/h with x/h as parameter bear a striking resemblance to 
the coefficient function shown in Fig. 3. Revised in accordance with the sill 
data now at hand, these diagrams are presented as Figs. 6 and 7, together 
with measured points for two representative values of x/h. The confocal rec- 
tilinear portions of the curves below h/w = 5 are identical with those previous- 
ly published, but the remainder will be found to differ appreciably. Similarity 
with the coefficient curve is maintained well into the sill range, however, ex- 
cept in the very neighborhood of the overfall limit. 

Although the foregoing experimental data were obtained with heads of rela- 
tively small magnitude and hence can be expected to apply only approximately 
to conditions at larger scale, it is believed that the degree of approximation 
will not be appreciably lower than the experimental precision evident in the 
figures. The Rehbock equation, for example, indicates that capillary and 
viscous effects change the discharge by about 1% at a head of 0.3 ft. Since 
this is representative of the heads that were investigated, there is no reason 
to expect scale effects of a higher order to be present in the results. 


CONCLUSIONS 


Through the combination of experimental measurements made at consider- 
ably different times and places, information is now at hand for the coefficient 
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of discharge and the form of the nappe profile as the relative height of a 
ventilated, sharp-crested weir at the end of a channel varies over its com- 
plete range from infinity to zero. In the range 0<h/w<6 the coefficient of 

weir discharge is given approximately by the formula 


h 
Ca = 0.61 + 0.08 rt 


In the range 0<w/h<0.06, the coefficient of discharge is given approximately 
by 


3/2 
C= 1.06(1 + :) 


In the intermediate zone there is a continuous transition between the two 
functions with a maximum at approximately h/w = 10; the latter might hence 
be regarded as the borderline between weir flow and sill flow. This value of 
h/w also marks the condition of maximum nappe deflection, the nappe for both 
weirs and sills gradually rising and thickening as this common value is ap- 
proached. Nappe profiles for design purposes may be determined either by 
interpolation of the profile curves or by direct reading of the generalized 
profile coordinates presented in the paper. 
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Discussion of 
“FREE OUTLET AND SELF-PRIMING ACTION OF CULVERTS” 


by Wen-Hsiung Li and Calvin C. Patterson 
(Proc. Paper 1009) 


WEN-HSIUNG Li,! A.M. ASCE, and CALVIN C. PATTERSON,2 J.M. 
ASCE—The writers are glad to receive discussion from investigators of 
three hydraulic laboratories. As their comments overlap to a certain extent, 
the authors would like to answer them together. 

The writers use mD for the mean piezometric head at the outlet in Eq. (1) 
and then define an effective mD by Eq. (4). As can be seen from the accom- 
panying sketch, and from Eq. (3) and (4), 


Effective mD = Actual mD + (4, 44) + (a, a 2) 


Note that the location of the projected hydraulic grade line corresponds to the 
effective mD, not the actual mD. Thus the values of effective mD presented 
in the paper are theoretically the same as those obtained by projecting the 
hydraulic grade line to the outlet section. This can also be seen from 


Pp 
Eq. (21) in Mr. French’s discussion by putting 4) Ge + z)dQ = actual mD, 
2 
Lv__ 
a .* Oo, and h, + 152g = A,. The small differences between 4, and 4, and 


between a and Op» are of no great importance in practice. Since A; and Ot; 


are unknown and will most likely be approximated in practice by Ao and a9 
respectively, the quantity effective mD is even more useful than the actual 
mD: 
2 


H + sL = effective 2g 


In presenting the data of effective mD, the writers call the quantity v/V gD 
a Froude Number. Mr. Blaisdell thinks that this quantity is not a Froude 
number “as ordinarily conceived” relating to free-surface flow. This is a 
matter of definition. It seems that the generally accepted definition of a 
Froude number is a dimensionless quantity in the form of v/YgD (or its vari- 
ations). It is an indication of the relative magnitudes of the inertia and gravi- 
tational forces(1) in a flow system, e.g. the flow from a closed container 
through an orifice.(2) Ina particular problem, the Froude number may some- 
times be interpreted in a particular way, such as the ratio of a characteristic 
velocity to the speed of propagation of an infinitesimal surface-wave, or as a 
measure of the ratio of the velocity head to the depth of flow or to the depth of 


1. Associate Prof. of Civ. Eng., The Johns Hopkins Univ., Baltimore, Md. 
2. Engr., Anderson and Nichols, Consulting Engrs., Baltimore, Md. 
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the conduit. Such interpretation may be helpful in particular cases but, it 
seems to the writers, should not be considered as a definition. The Froude 
number is significant in determining the value of mD because the pressure 
distribution at the outlet section depends on the relative magnitudes of inertia 
and gravitational forces acting on the jet. The writers use v/VgD instead of 
one of its variations to take advantage of the fact that, at the same v/ VeD, the 
values of m for unsupported jets are about the same for square and circular 
pipes. 

Considering the nature of the measurements, the differences between 
Mr. Flack’s data and the writers’ data of m for horizontal culverts with no 
outlet-apron are not unexpected. It is believed that the data presented for 
various culvert slopes and outlet conditions are consistent enough for practi- 
cal purposes. 

The writers appreciate the comments on the self-priming action of cul- 
verts. The work of Mr. French in this respect is particularly interesting to 
the writers. They are glad that Mr. French adds Type A to the three types of 
priming mentioned by the writers who, in writing the paper, considered 
Type A as the limiting case of Type 1. It seems that Mr. French’s Type B 
of priming is the transient phenomenon of what the writers call Type 2 prim- 
ing. In the paper, the writers suggest a way of presenting the data on priming 
for practical use. As there is no exact knowledge about turbulent non-uniform 
flow, air entrainment and vortex formation (as of 1957), there is no exact way 
of transferring model data to a prototype. A practicable, though approximate, 
method has to be used. It should be recalled that, in using any empirical in- 
formation, it is a sound engineering practice to allow a margin of safety. The 
data presented in the paper, being experimental results, should be used in the 
same spirit. Friction is a major influence in Type 1 and Type 2 priming. In 


Type 3 priming where air entrainment is most important and friction only 
1/3 
plays a secondary part, the action depends mainly on D/H. The use of a 


is not incorrect though not absolutely necessary. 
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Discussion of 
“THE APPLICATION OF SEDIMENT-TRANSPORT MECHANICS 
TO STABLE-CHANNEL DESIGN 


by Emmett M. Laursen 
(Proc. Paper 1034) 


EMMETT M. LAURSEN.!~—A discussion such as that of Mr. Shulits can 
only gladden the heart of a writer, for it amplifies and reinforces the original 
contentions set forth. The ‘formula dilemma,” as Mr. Shulits puts it, is such 
that at the present time the design of stable channels cannot proceed with any 
substantial degree of confidence on the basis of sediment-transport mechanics 
alone—experience is absolutely essential. Sediment-transport mechanics 
can serve as a guide and the many formulas can serve as scaling relation- 
ships, however, and in that way supplement experience. When the formulas 
are used to compare different channels the gross quantitative disagreements 
largely disappear and the similarities are emphasized. 

In a personal communication Mr. J. W. Delleur has pointed out that power 
of n in the Brown-Einstein equation is incorrect. In Table I the power should 
be 3 and in Table II the power should be 6/5. 


1. Research Engr., Iowa Inst. of Hydr., Research, Iowa City, Iowa. 
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Discussion of 
“MEASURING STREAMFLOW UNDER ICE CONDITIONS” 


by A. M. Moore 
(Proc. Paper 1162) 


CHARLES E. BEHLKE,* A.M. ASCE.—As more and more of the cheaper 
water power sites in the United States are developed, the day approaches 
when American industry must look to the great power potential of northern 
Canada and the Territory of Alaska. In these areas of great water power 
potential the winters may last as long as eight or nine months each year. It 

is in the North that the methods of evaluating ice-affected records of stream 
flow, so well described by Mr. Moore, have greatly magnified significance. 

It is quite common in Arctic and subarctic regions for the winter stream- 
flow to constitute a relatively minor amount of the total annual flow. How- 
ever, since stream beds may drop quite rapidly, the winter flow can be quite 
important to a hydro development because of the frequent’lack of potential 
reservoir volume. It is therefore quite important that streamflow records be 
obtained for the winter months. With a very few exceptions, the methods 
which Mr. Moore has so adequately described of evaluating streamflow from 
stage records would be quite applicable to these northern areas. 

Discharge measurements by themselves are quite tedious and difficult to 
make in cold weather especially if the thickness of the ice cover on the 
stream is very great. It is quite common for ice thicknesses on northern 
streams to range from Zero to six feet. This presents a considerable prob- 
lem before the current meter can be placed in the stream. The efficiency of 
a man gets quite low as the temperature drops, so a man taking a discharge 
measurement in the very cold weather of the winter may take several days to 
obtain a measurement which would take him only a few hours in the warmer 
weather. To add to the difficulty, the daylight hours are much shorter in 
these northern areas than they are further south, so a single discharge mea- 
surement may take two men two or three days to complete if the weather is 
not too cold or otherwise inclement. If the weather becomes colder than 
about minus 40°, very little is lost by simply stopping work and waiting for 
warmer weather. This means that men must wait at a location sometimes for 
long periods of time before they are able to make the discharge measurement. 
It is not uncommon in Alaska and northern Canada for the temperature to drop 
to as low as minus 60°, and colder, which means that during a normal winter 
there is a good deal of time during which the weather is too cold to work out- 
side efficiently. 

Since discharge measurements during the winter in northern climates are 
quite expensive, streams may be gaged only a few times during these months. 
The discharge is then known only at the time of the measurement and must be 
estimated for other times. A record of the stage, even though it is effected 
by ice cover, anchor ice, or ice jams, would be quite valuable for a deter- 
mination of the discharge during winter months as has been shown by 
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Mr. Moore’s paper. The obtaining of such a record in a cold climate is, at 
present, practically impossible because of the limitations of the methods 
used for obtaining these records. The standard type of gage well as well as 
the intake pipes leading from the river to the well freeze quite early in the 
fall and remain frozen quite late in the spring and even into the summer. The 
writer has had the experience of being a member of a party which thawed the 
intake pipes to a gage house on the Tanana River, one of Alaska’s major 
rivers, after the first of July. This freezing can affect spring and summer 
records which are very important. In Alaska, since the frost has usually 
penetrated deepest into the ground sometime in May, intake pipes thawed af- 
ter the spring runoff begins may refreeze again very quickly. This points up 
the very serious fact that considerable changes must be made in the matter 
of obtaining stage relationships for these northern streams before good re- 
cords can be obtained. 

The writer would warn the inexperienced engineer against attempting to 
correlate stream discharges with air temperatures at low temperatures. 
Mr. Moore has pointed out that this can be done at moderate temperatures 
below freezing. The writer has seen engineers very carefully adjust their 
estimates of stream discharge with changes in temperature where the tem- 
perature ranges were in the order of from minus 20° to minus 65° with the 
cover on the ground perhaps averaging 3 feet of very light, well insulating 
snow. Since winter stream flow in northern climates is derived entirely 
from groundwater, it is quite improbable that the particular day’s discharge 
from a stream could be affected by the temperature for that day, but it would 
be more a function of a mean temperature extending over a much longer 
period of time as well as the depth and density of snowcover on the basin. 

Some engineers make a practice of installing thermograph thermocouples 
in the gage well instead of in the stream. This can make a considerable dif- 
ference in the temperature which would be recorded for the stream at a par- 
ticular time. The writer has measured differences between gage well and 
stream temperatures of upwards from fifteen degrees F. in the summer. In 
the winter after the gage freezes, the temperature for the stream would 
erroneously be recorded as below freezing by perhaps a considerable amount. 
This is obviously not a true record of stream temperatures, and the use of 


such records should be reviewed critically by the engineer having recourse 
to them. 


STEPONAS KOLUPAILA.*—This paper by Mr. Moore is a valuable contri- 
bution to the complicated problem of determination of streamflow under ice 
conditions. The phenomen of siphon action is of particular interest. The 
writer had occasions to observe similar strange effects in a cross section 
jammed with frazil ice: the water flow, partly in the form of a slurry, was 
concentrated in a few separate channels under extremely different hydraulic 
conditions. There is no substitute for intelligence and experience in compu- 
tations of winter flow, more so than in any other field. 

The writer was very interested in the method of adjustment of daily dis- 
charges applied by the author. Obviously, the author uses the Stout method: 
he corrects the water stages to the “effective heights,” as this is explained, 
e.g., in the “Stream Flow” by Grover and Harrington, p. 285. Another method 
of correction of discharges seems to be of greater flexibility and accuracy: 


* Prof. of Civ. Eng., Univ. of Notre Dame, Notre Dame, Ind. 
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the ratio is determined between the actually measured water discharge Q' 
and that read on the discharge curve (available for open channel conditions) 

Q. This ratio k = Q' /Q varies through the winter period due to the changes in 
thickness of ice cover, presence or disappearance of frazil ice, depth of water 
and, particularly, the smoothing of subsurface roughness by flowing water. 
The obtained ratio is to be interpolated between actual measurements, ad- 
justed at the beginning and the end of ice cover, when direct measurements 
are too dangerous or impossible; meteorological data can be useful in tracing 
the diagram of k. Daily discharges are to be computed multiplying the value 
read from discharge curve and the ratio k from diagram. Fig. 1 illustrates 
this method. 

This “discharge correction method” in contrast to the “gage correction 
method” was initiated by the writer in 1928 and since that time has been in 
use in many European countries. Several detailed procedures were recom- 
mended for deficient data, especially for previous years without any direct 
winter measurements. Thus, the writer was able to compute the runoff of the 
Nemunas River in Lithuania for 1812 - 1943, for complete 132 years, one of 
the longest series of existing reliable runoff data. Further developments of 
this method might be mentioned. Mr. A. Ogievskil derived a “hydrometeoro- 
logical method” for corrections k: air temperature is used; accumulation of 
the sub-freezing temperctures aids in construction and detailization of the 
diagram. Recently, Mr. H. A. Klein obtained an approximate empirical formu- 
la for the variation of the k ratio in Elbe River: k = 0.308 + 0.02 T, where T 
is the number of days since ice cover was established. 

This method was presented by the writer at several international confer- 
ences: Leningrad 1928, Tallinn 1928, Seville 1929, Barcelona 1929, 
Helsingfors 1936, Edinbourgh 1936, Lubeck 1938. In this country Mr. N. C. 
Grover, chief engineer in the Water Resources Division, U. S. Geological 
Survey, tried in 1936 to introduce this, as he called it, “Lithuanian method” 
into practice of his Division. Several valuable reports were submitted on this 
subject, as by M. C. Boyer and B. J. Peterson, G. E. Harbeck, C. C. McDonald 
and others. Unfortunately, these interesting reports were not published in the 
Water-Supply Papers. The writer is indebted to Mr. McDonald for a copy of 
his dissertation “Determination of stream flow for periods of ice effect,” 
1943, which is of very great value and still retains its importance. 

The writer disagrees with the statement of Mr. Moore that engineering 
literature contains few references concerning the effect of ice on stream flow 
measurements. Starting with the pioneer work, done by eminent American 
hydrologists H. K. Barrows, R. E. Horton, and W. G. Hoyt, there are many 
papers in many languages. Some of the most interesting are presented in the 
subsequent reference list. 

The opinion of the author concerning the accuracy of winter records is 
completely true. His suggestion to reexamine all previous winter flow compu- 
tations before 1950 and to improve the annual runoff data is of great signifi- 
cance and should be materialized. The publication of reports on winter flow 
would contribute much to the progress of establishment of the runoff balance. 
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Discussion of 
“TECHNICAL PROBLEMS OF FLOOD INSURANCE” 


by H. A. Foster 
(Proc. Paper 1165) 


G. N. ALEXANDER, * A.M. ASCE.—Mr. Foster’s latest paper concerning 
the application of statistics to flood estimation entitled “Technical Problems 
on Flood Insurance” is, as he states, based largely on his preceding paper on 
“Flood Insurance.”(1) His first article on statistics for engineers, published 
by the ASCE some thirty years ago(2) still remains a classic in this field. 
As the later papers give insufficient technical background to the author’s 
findings concerning flood probabilities it is not possible in this discussion to 
enter into the technical detail that characterised his first article. 

The two basic questions concerning flood probabilities are: 


Question (A) What is the “best estimate” of the probability or return 
period of a flood greater than a specified magnitude? 
To answer this question decisions are required on— 
(1) The type of distribution 
(2) The method of parameter estimation. 


Question (B) What is the “precision” or reliability of this estimate? 


The earlier papers, such as those of Foster(2) and Goodrich(3) dealt only 
with the first step of Question A, viz. Distribution Type. There is little dis- 
cussion as yet on the second step in engineering journals, but it is treated in 
statistical text books. 

The importance of Question B was recognised some thirty years ago by 


the actuary Arne Fisher, in his discussion on Goodrich’s paper(3) when he 
stated— 


“the paper deals only with a subject of graduation, and leaves almost 
untouched the far more subtle and difficult aspects of the probable or 
presumptive values of future observations.” 


As statistics is applied to flood data primarily to enable inferences to be 
made concerning the magnitude of future floods, this problem is paramount. 
This question of precision is closely related to that of sampling errors, which 
as the author states is now being recognised by the engineering hydrologist, 
although as yet no particular procedures are generally accepted. 

Given sufficient data and using recognised statistical techniques, Ques- 
tions A & B above can be answered for recurrence intervals not greatly ex- 
ceeding the length of record; the difficulties arise when extrapolation is 
necessary. For an assumed distribution, the precision of an estimate de- 
creases approximately as the logarithm of the recurrence interval and in- 
creases as the square root of the length of record. Consequently extrapola- 
tion of short records can give quite misleading results, but with long records 


* Engr.-in-Charge of Water Resources Section, State Rivers and Water Sup- 
ply Comm., Melbourne, Australia. 
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some extrapolation is possible without undue loss of accuracy. 

Failure to recognise the magnitude of the sampling errors invariably 
associated with short records was largely responsible for discrediting the 
use of probability methods in relation to floods. The writer believes that 
statistical techniques properly applied must ultimately form the basis of es- 
timating flood probabilities which in turn are essential to actuarial estimates 
of premiums for flood insurance. 

We next examine more closely, but still very generally, the means of re- 
ducing the arbitrariness in current procedures in seeking to answer Ques- 
tions A and B. 


Question (A) Methods of Curve-Fitting 
As the author states “the various proposed methods of curve-fitting 
generally produce reasonably good results with a given flood record when 


used to express probabilities or recurrence intervals within the length of 
record.” 


Step (1) Type of Distribution: 

Attempts have been made by some authors to justify certain types of dis- 
tribution for flood estimation on “a priori” grounds. Thus, floods are re- 
garded by Gumbel(4) and others as “extreme values” and hence they favour 
the Fisher-Tippett type of distribution. Others contend that the log-normal 
distribution is basically more appropriate, as floods may be regarded as the 
product of a number of independent causes. Arguments of this sort carry 
little weight, as they are seldom backed by study of causes of floods and 
their hydrographs; also the conventional use of the Chi-squared test for 
“goodness of fit” is of doubtful utility. 

Chow has shown(5) that for a specific value of the variance, the asymptotic 
extreme value distribution and log-normal distribution when standardized are 
practically identical. From unpublished analyses carried out in the writer’s 
office, it appears that the log-normal can also be used in lieu of the finite ex- 
treme value distribution based on a normal parent by selection of the appro- 
priate parent size. This suggests that the log-normal distribution is suffi- 
ciently “flexible” for most cases. The statistical characteristics of the 
log-normal are becoming better known and its conservatism in extrapolation 
renders it particularly suitable for estimating design floods for spillways. 


Step (2) Estimation of Parameters: 

As long as there are insufficient grounds for differentiating between par- 
ticular distributions and methods of parameter estimation, there is some 
justification for the time-saving process of fitting curves “by eye.” However, 
if extrapolated values are quoted, it is desirable to specify the type of distri- 
bution used and the method of fitting, for the different procedures give at 
times very different answers. 

Having selected the distribution, a decision is now required on the most 
appropriate method of estimating parameters. Two methods are commonly 
used by statisticians viz. the “method of moments” and method of “maximum 
likelihood.” For certain theoretical reasons the maximum likelihood method 
is preferred by most statisticians, at least for fitting within the range of ob- 
servations. However, if curves have to be extrapolated the method of mo- 
ments automatically weights the extreme values, and there may on this ac- 
count be some justification for its use—but this is a contentious technical 
problem outside the scope of this discussion. Even when fitting a log-normal 
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distribution, estimates from extrapolated curves will differ significantly de- 
pending on whether the fit is to the original observations, or—as appears 
more desirable—to their logarithms. 


Question (B) Precision of Estimates 
The author refers to the question of the reliability or precision of esti- 
mates towards the end of the paper stating “the question of the reliability of 
probability estimates of future floods is of much greater significance in a 
programme of flood insurance than in justifying a flood protection project.” 
Any estimates are of much greater value if some indication is given of their 
reliability and this is particularly true of flood estimates. 
This question of the precision of estimates as a function of sample size is 
a basic problem of statistical inference. It may be said to date from “Stu- 
dent’s” work on the t-test published nearly 50 years ago. However, there is 
still controversy concerning the more technical problems of statistical in- 
ference such as the relation between the Neyman-Pearson concept of “con- 
fidence limits” and Fisher’s “fiducial limits?”(6) In addition to these two 
methodological developments, “tolerance limits” may also be obtained over 
the range of observations. Hydrologists seeking to apply these procedures 
concerning probability limits naturally have some doubts in trying to decide 
between them—apart from any mathematical difficulties that may arise. 
Although the practical differences in the results using alternative methods 
may be small within the range of observations, the validity of inferences over 
the extrapolated estimates is questionable, particularly when the fitted curve 
depends on the distribution type assumed. For example, when the log-normal 
distribution and the Pearson Type III were both fitted to 50 years of annual 
maxima of monthly discharges of the River Murray at Jingellic, it was found 
that the log-normal estimate for 100 year return period was very close to the 
upper 2 1/2% risk curve (95% confidence limits) for the Pearson Type II. 
Such differences, although disturbing to those unfamiliar with the character- 
istics of statistical methods, are believed to be no greater than differences 
associated with other more popular methods of floods estimation for long 
recurrence intervals such as the “maximum possible precipitation” technique. 
The rather misleading title to this technique implies that no particular 
probability is associated with the estimate obtained, nor is its precision in- 
dicated. The method does show how, given certain meteorological conditions, 
floods can occur whose recurrence intervals are so great that probability 
methods which ignore the effect of sampling errors need re-interpretation. 
It is believed that this meteorological method and that of statistical inference 
from flood data will be reconciled when the probabilities associated with such 
steps as storm maximization and transposition can be ascertained. 


CONCLUSION 


It is clear that there is still much to be done both statistically and hydro- 
logically before the estimation of future floods by statistical inference is 
placed on a firm foundation. With the recent enactment of flood insurance 
legislation in U.S.A. and the consequent need for actuarial determination of 
premiums, it may be expected that studies of this sort will be accelerated. 
Mr. Foster has made notable contributions in this field, and it is hoped that 
this paper will stimulate others to publish their findings. 
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STEPONAS KOLUPAILA.*—In his discussion of the determination of flood 
probability, the author states: “After studying and comparing various 
methods of curve-fitting that have been proposed for hydrological studies, the 
author has not been able to reach any definite conclusions regarding which 
one is the best for use in flood-insurance studies.” As it is known, Mr. 
Foster offered his method in 1924 and prepared his tables, based on two 
characteristic values—coefficient of variation and coefficient of skew. Fos- 
ter’s method was quite popular and widely used here and overseas. After 
some unfortunate failures occurred in extrapolation of probability curve this 
method was almost entirely eliminated. The writer feels that the cause of 
such a failure was the improper application and lack of experience. The co- 
efficient of variation can be easily established and checked, although it may 
be based on short period of available data. The coefficient of skew is more 
intricate: the usual period of 20 to 30 years is definitely not sufficient for 
establishment of any reasonable value. Sometimes the addition of just one 
year causes a great change of this coefficient. Therefore it is more sound to 
presume a value by analogy with some more reliable long range data. The 
adjustment suggested by A. Hazen for inadequate length of records is not 
satisfactory. Different rules were offered for improvement of the coefficient 
of skew. Russian hydrologists made some significant progress: they intro- 
duced different empirical relationships and an additional “guarantee-factor.” 
The results obtained with these corrections are more satisfactory. The 
writer has a successful experience with Foster’s method and therefore hopes 
that this excellent American method will be restored at home in near future. 
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Discussion of 
“FREQUENCY ANALYSIS OF STREAMFLOW DATA” 


by David K. Todd 
(Proc. Paper 1166) 


CORRECTIONS—Through an oversight, Figs. 10 through 16 were not pub- 
lished with this paper when it was published in February, 1957. These 
figures, therefore, are reproduced herewith. 


M. A. BENSON.!—The author is to be commended for attempting to bring 
before engineers a summary of information on the application of statistical 
analysis to streamflow data. The subject is not entirely new to engineers. 
For a comprehensive discussion of flood-frequency practices and a complete 
bibliography on the subject the writer would like to call attention to the recent 
article of Chow2 on this subject. 

Mr. Todd refers to available data in the annual Water-Supply Papers of the 
U. S. Geological Survey. Mention should also be made of the compilations of 
records of surface waters of the United States through September 1950, now 
being published in a series of volumes, each of which covers a major drain- 
age basin of the United States. Volumes already published are in Water- 
Supply Papers 1301, 1311, 1313, 1316, and 1317 respectively. These include 
listings, in readily available form, of annual peak discharges. 

The author states that, “For temporary or small structures where designs 
are based upon flood frequencies about 10 years or less, the basic flood 
method (partial-duration series) should be used;---.” The author repeats 
here the recommendation of the Sub-committee of the A.S.C.E. Joint Division 
Committee on Floods. However, experience has shown that it is preferable in 
most cases to use the annual-flood method, even for small design frequencies. 
The annual-flood method is simpler, the problem of the independence of 
floods is not then involved, and it is always possible to convert from the 
annual-flood to the partial-duration series (see Langbein, reference (22) in 
original article). 

The author says “Most extrapolation by eye of frequencies is arbitrary, 
hence it does not provide a sound basis for estimating frequencies of high 
floods. To overcome this difficulty analytical techniques have been developed 
for fitting a cumulative frequency curve to specified types of frequency dis- 
tributions.” He says further that “With the aid of statistical parameters de- 
termined from the streamflow data and special plotting papers, it is possible 
to estimate say the 10,000-year flood from a record of 20 years duration. 
That such practices are dangerous has been well emphasized by the ASCE 
Subcommittee on Review of Flood Frequency Methods.” 


1. Hydraulic Engr., U. S. G. S., U. S. Dept. of the Interior, Washington 25, 
D. C. 

2. Ven Te Chow, Hydrologic Studies of Floods in the United States, Symposia 
Darcy, Dijon 1956-Tome III, Union Geodesique et Geophysique 
Internationale, p. 134-170. 
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FIG. 12 DAILY, WEEKLY AND ANNUAL-FLOW DURATION 
CURVES, BIG WALNUT CREEK AT REES, OHIO. 
(AFTER 10) 
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FIG. 14 OBSERVED AND FITTED FREQUENCY DIS- 
TRIBUTIONS OF MEAN DAILY FLOWS, 
CUCAMONGA CREEK AT UPLAND, CALIF. 

(AFTER 27) 


1200 | 

200 
-1.0 1.0 20 30 


DISCUSSION 


O 


> 
O 
uJ 
uJ 


LOG DISCHARGE 


FIG. 15 FREQUENCY DISTRIBUTION OF MEAN 
DAILY FLOWS, FISH CREEK NEAR 
DUARTE, CALIFORNIA. (AFTER 27) 
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FIG. 16 FREQUENCY DISTRIBUTION OF MEAN 
DAILY FLOWS, SALMON RIVER AT 
SOMESBAR, CALIFORNIA. (AFTER 27) 
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The writer questions the implied superiority of fitted mathematical dis- 
tributions over graphical techniques. It has not yet been possible to demon- 
strate a priori that flood peaks should conform to any mathematical or sta- 
tistical distribution. Hence, it is not only dangerous but meaningless to 
estimate a 10,000-year flood from a 20-year record by means of a fitted 
curve. Pertinent statements on the general subject of graphical versus 
mathematically-fitted curves have been made by Ezekiel:3 


“When there is some logical basis for the selection of a particular 
equation, the equation and the corresponding curve may provide a definite 
logical measurement of the nature of the relationship. When no such logi- 
cal basis can be developed, a curve fitted by a definite equation yields on- 
ly an empirical statement of the relationship and may fail to show the true 
relation. In such cases a curve fitted freehand by graphic methods, and 
conforming to logical limitations on its shape, may be even more valuable 
as a description of the facts of the relationship than a definite equation and 
corresponding curve selected empirically. 

“In any event, estimates of the probable value of the dependent variable 
cannot be made with any degree of accuracy for values of the independent 
variable beyond the limits of the cases observed; and can be made most 
accurately only within the range where a considerable number of observa- 
tions is available. It may be possible to extrapolate the curve if its equa- 
tion is based on a logical analysis of the relation as well as on the cases 
observed; but in that case the logical analysis, and not the statistical ex- 
amination, must bear the responsibility for the validity of the procedure.” 


See also Thomas,4 who says, “Indeed insofar as extrapolation made be- 
yond the range of the recorded floods is concerned, the most elaborate analy- 
tical procedure of curve-fitting gives results that are no more reliable than 
those obtained by a simple estimation by eye of the flood frequency curve on 
any kind of probability paper.” 

The author refers briefly to the method of regional flood-frequency analy- 
sis which can be employed to estimate flood frequencies in ungaged areas. 
His figure 6 illustrates such methods, developed for the Atlanta, Ga. area, but 
he omits mention of the many statewide regional-frequency studies in which 
such methods have been described and developed for a large part of the 
country. 

The author discusses the fact that in some regions floods may be caused 
by either rain or snow melt and, because ‘two different populations are inter- 
mingled----- it follows, therefore, that not one but two straight lines of differ- 
ing slope should be used to fit such flood frequency curves.” Because the 
populations are in truth intermingled, the use of two straight lines does not 
follow as a solution. There is in general no particular flood magnitude below 
which all floods are caused by one factor and above which they are caused by 
another. In most parts of the country floods have multiple causes which are 
not distinct but combine or merge one into the other. The situation has been 
well analyzed by Hazen,° who many years ago analyzed the effect on 


- Mordecai Ezekiel, Methods of Correlation Analysis, 2nd Edition, 1941, 
John Wiley & Sons, N. Y., p. 127. 
4. Harold A. Thomas, Jr., Frequency of Minor Floods, Journal, Boston 
Society of Civil Engineers, October 1948, p. 428. 
. Allen Hazen, Flood Flows, John Wiley & Sons, N. Y., 1930, p. 102. 
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frequency curves of the multiple causes of floods, but said, “It is true of 
course that these causes grade into each other and that no lines of demarca- 
tion can be drawn.” 

The discussions of the frequency analysis of low flows and of all flows 
could be expanded considerably, if a survey of these subjects is to be at all 
complete. 
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Discussion of 
“BUTTERFLY VALVE FLOW CHARACTERISTICS” 


M. B. McPherson, H. S. Strausser and J. C. Williams, Jr. 
(Proc. Paper 1167) 


TURGUT SARPKAYA,! J.M. ASCE.—Valuable information adding to the 
gradually increasing knowledge of the performance of butterfly valves is con- 
tributed by the authors. The paper constitutes a well-arranged review of 
previous studies in addition to the newly presented basic analytical relation- 
ships, their interpretations and sensible applications to actual flow condi- 
tions. However, the authors have not dealt with the following features of the 
problem: (a) Variation of contraction coefficient at each end of the valve 
plate with the angle of valve opening and with the angle of complete closure; 
(b) Stability of flow for various valve openings; and (c) An explanation for the 
large scatter of experimental data for small values of o& (@<20°) or for near- 
full-open positions of valve plate. 

For some time, this writer has been interested in tne application of mathe- 
matical methods to practical problems of fluid flow, and has contended that 
many existing semi-empirical relationships would find wider application if 
all parts of the relationships were, wherever possible, freed from empirical 
values and experimental coefficients. Having once obtained the pure mathe- 
matical results, engineers can then compare them with the experimental re- 
sults. Since the writer feels that the results obtained by the authors are sig- 
nificant, it seems worthwhile to apply the methods of hydrodynamics to 
butterfly valve flow in order to have a better understanding of the variations 
of various parameters involved in the phenomenon. From this viewpoint, the 
omission of any discussion of the rather surprising analogy which permits the 
use of hydrodynamic theory for certain problems of viscous motion is deemed 
to be regrettable. Having this object in mind the theoretical solution of the 
present problem is given in this discussion. Then the experimental results 
are compared with the theoretical ones and answers provided for the points 
stated above. 

It is well known that many plane irrotational flow patterns which include 
the formation or deflection of free jets can be analyzed completely by the 
method known as the Helmholtz(1)-Kirchhoff(2) theory of free streamlines. 
The results obtained may be significant not only for the corresponding two- 
dimensional flows of real incompressible fluids but, in some instances, for 
their three-dimensional counterparts as well. The analyses consists of the 
definition of successive conformal transformations involving a hodograph, or 
velocity plane, and the application of Schwartz(3)-Christoffel(4) transforma- 
tion. The writer would like to emphasize the fact that from the mathematical 
viewpoint, the direct calculation of free-streamline flows is restricted to two- 
dimensional irrotational flows of incompressible fluids which are free from 
gravitational effects and for which the solid boundaries are straight. As flows 
of real fluids seldom satisfy fully these stringent restrictions, the practical 
value of the results is open to question. The significance of the various 


1. Nazilli, Turkey. 
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restrictions can be shown only by comparing the calculated results with those 
obtained from observations of real flows. For the present problem since the 
shear is not a primary factor, (the effect of shear will be discussed later) 
for the flow accelerates rapidly, the differences between comparable ideal 
and real flows are small. The importance of gravitational effects for the 
coefficient of contraction is already known to be definitely secondary. Being 
fully aware of the restrictions and limitations of the applicability of this 
method, the solution of butterfly valve problem proceeds as follows: If the 
coefficients of contraction at points B and H are Cc and Ce2 respectively as 
is shown in the physical plane (Fig. 1), then the main objective of the present 
investigation becomes the determination of the relationships: 


The existing relationships are: 


(3) 


(4a),(4b) 


case: 


(2) 
since 
m= = 
C..+ Cos (5) 
sing 
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in which n-and m-represent the thickness of the jets at B-and H-respective- 
ly; 7 is the width of channel; B the angle of closure; @ the variable angle; and 
L the length of valve plate. The velocities of both jets are equal to Vj and the 
velocity far upstream from the valve plate is Vp. 

If z-and w-are used for the complex variable and the complex potential, 
then the variable £ defined as (dw/dz) is given by 


(6) 
S=-u+iv=—qge ie 


in which u-and v-are the velocity components in the x-and y-directions, q is 
the magnitude of the velocity vector, and @ is the angle of its inclination. 
If the variable 2, 


(7) 


is introduced, the flow region in the physical plane (Fig. 1) can be mapped on- 
to a corresponding region in the Q -plane, as indicated in Fig. 2. Either ¢ -or 
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Q-must be related to w by means of known methods through still another com- 
plex variable designated as t. This relationship completes the rather com- 
plicated functional analysis because, as has already been seen, either ¢ -or 
Q-can be expressed in terms of (dw/dz). 

Any polygonal boundary in one complex plane can be transformed into the 
entire real axis of another through application of the Schwartz-Christoffel 
theorem, the interior of the polygon being transformed into either the upper 
or the lower half of the second plane. Thus, the outline in the 2-plane is 
directly transformable into the real axis of the t-plane. Therefore, if we open 
out the Q2-plane and lay it along the real axis of the t-plane as shown in Fig. 3, 
the interior of the 2-plane transforms into the lower half of the t-plane. The 
transforming function can be written as follows: 


Ht 
ace 
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dt 


(tt, 


t. Plane 
FIG.3 


As many factors (t - ty) are introduced into the transformation as there are 
vertices of the polygonal boundary in the 2-plane. The value tp is the locus 
of the transformed point on the real axis in the t-plane, and ¥, is the exterior 
angle of the polygon at that point in the t-plane. Three of the tn values can be 
assigned arbitrarily but with due regard to the geometrical characteristics of 
the Q-plane, and the remainder must be given parametric values to be evalu- 
ated in terms of given quantities, along with M-and N, from the integrated 
function. 


The desired mapping function between Q2-and t-becomes: 


(t+a) dt 
=™M 


The transformation equation between the w-plane and the t-plane is ob- 
tained directly from the complex potentials for sources and sinks, point A be- 
ing considered as a source of strength (27)V,) and points B and H as sinks of 
strengths (2nVj) and (2mVj) respectively: 


{nin(t-1) +m ln (t+1) (10) 


wal + G (G+ G,)Ln(t+a)} (10a) 


and w-is connected to z by the relation, 


(11) 
+/e 


The evaluation of this and the other foregoing integrals are shown in the 
appendix of this discussion paper in order not to lose sight of the practical 
nature of the problem undcr consideration. 

The results or the functional relationships obtained from these equations 


1348-34 

-f -l +1 d 

F H A B OD E, 

| 

or 


ASCE 


+ 

ean 


DISCUSSION 


4-4-4 

sen + 
enone ees 
ttt +1 

tr 


= 
wes 
age: 


are shown in Fig. 4. The significance of this figure is obvious. The contrac- 
tion coefficient C-j varies greatly with a, but very little with the angle of 
closure 8 , whereas the contraction coefficient Ceg varies appreciably with 6, 
but very little witha, for the values of a greater than approximately 15°. 

But for small values of a, Cc2 varies extremely fast and approaches unity. 
The values of the corresponding angles (corresponding to@<10°) used in 
evaluating the elliptic functions (see appendix) approach their limit so closely 
that difficulties are encountered because of limitations in the tables available 
to the writer.(5,6,7,8,9,10) For those portions of the curves which could not 
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be computed readily the trend is indicated by broken lines. 
In accordance with the authors’ equations (1) and (7a), the discharge coef- 
ficient Co for free discharge can be written as follows: 


Cg values computed from Eq. (12) for B = 75° using the contraction coeffi- 
cients obtained in Fig. 4, are plotted in Fig. 5, and the solid-line is obtained. 
In the same figure the experimental data obtained by the authors for 4"'CDC 
valve and 6'"'CDC valve are also plotted. Though, the complete closure angle 
of the latter valve differs 2.5° from the former, another theoretical curve has 
not been drawn for 8 = 77.5° in order not to make the figure illegible. The 
agreement between the theoretical and experimental results is fairly good, 
the theoretical discharge coefficient being somewhat smaller. 

A similar computation has been made for 8 = 90°, and the resultant the- 
oretical curve with the experimental points (Genissiat and Dow data, as re- 
ported by the authors) are shown in Fig. 6. The agreement between the the- 
oretical and experimental results is excellent. 

So far we have discussed the effects of geometric variables on each coef- 
ficient of contraction and on the discharge coefficient. It is seen that the 
coefficients of contraction are functions of a-and B-only. Under the heading, 
“General Analysis for Enclosed Flow,” authors state that: “In these tests, 
the magnitude of y (see authors’ Fig. 1) was not noticeably changed by varying 
the discharge.” Authors are perfectly correct in their observations, because, 
theoretically the value of y (n-and m-in Fig. 1 of this discussion paper) does 
not depend on the magnitude of discharge. 

The second point we have proposed to discuss was the stability of flow for 
various valve openings and in connection with this the explanation for large 
scatter of experimental data for very small values of a@-or for near-full-open 
positions of valve plate. The stability of flow depends on the stability of the 
stagnation point on the valve plate. As a@-becomes smaller, the stagnation 
point or rather specifically speaking the stagnation line approaches the lower 
edge of the valve plate and theoretically it reaches the edge when a-is zero. 
But, practically the stagnation point or line reaches the sharp end of the plate 
when a@- becomes about 15°. To illustrate this, a curve has been drawn com- 
puting the location of the stagnation line for 8 = 90° and for various values of 
a@, as shown in Fig. 7. The dimensionless parameter A -becomes practically 
zero for a = 20°. 

In the case of axially symmetrical flow, the locus of stagnation points is 
curvilinear. Therefore, the stagnation points in the vicinity of valve axis are 
very close to the edge of valve plate, whereas the stagnation points in the 
vicinity of the diameter perpendicular to the valve axis are relatively far 
from the edge. The writer believes that in the case of axially symmetrical 
flow for a critical value of a-the stagnation points relatively near the plate 
edge jump abruptly to the edge of the plate changing the discharge coefficient. 
In the two-dimensional case the entire stagnation line jumps abruptly to the 
edge of the plate for small values of @. For a critical value of a, however, 
the stagnation line is unstable and jump to the edge and then back alternately 
resulting in unstable flow. 

It would be very interesting to compare the discharge and contraction 
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coefficients obtained from two separate tests as follows: (a) First setting 

a = 250 and decreasing it a few degrees each time, (b) First setting a@ = 5° 
and increasing it a few degrees each time. Each test might well result ina 
separate curve in the range of critical values of @. The results obtained 
from such an experiment should prove useful in obtaining a better understand- 
ing of the stability of flow. 

It seems to be evident from the foregoing explanation that the scatter of the 
experimental data for small values of a-is largely due to the rapid and abrupt 
variations in the location of stagnation line in that region. Moreover, the 
geometry of valve plate also has an important effect on the variations of afore- 
mentioned coefficients. Theoretically, the plate is a thicknessless or one 
dimensional surface, but practically it has a thickness. Therefore, as noted 
by the authors under the heading: “Comparison of Actual and Analytical Flow 
Coefficients,” “For a blade with a relatively wide tip, the point of separation 
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will be at one or the other edge, depending upon the relative blade angle. The 4 
sudden, slight shift of Co trends for some valves is due to this transfer of the 
separation point.” 


The writer would like finally to comment on the discharge and on the con- 
traction coefficients for submerged flows. It is evident from the experimental 
data given by the authors that the discharge coefficients for submerged flows 
are larger than the ones for free discharge as amply explained by the authors. 
They have observed that as stated under the heading: “Free discharge and 
submerged flow characteristics,” “The free-discharge and submerged char- 
acteristics appear to be related, but are distinctly different at near-open 
positions. In submerged flow, the fluid fills or nearly fills the void between 
the two jets up to the blade, resulting in an expansion of the given fluid within 
the same fluid. In free-discharge of a liquid, the void is normally filled with 
air. For this reason, it would be logical to assume that the flow characteris- 
tics might not be identical even though the installations are quite similar.” 

The authors, therefore, assume that the action of turbulent shear in 
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causing the submerged jets to expand has some influence upon the pressure 
field that controls the rate of efflux. This writer was inclined to accept this 
explanation if it were not due to another contradictory observation and state- 
ment made by Mr. Curtis(11) and his associates on the “Coefficient of con- 
traction for a submerged jet.” Mr. Curtis states that “The minimum dimen- 
sion of the submerged-jet profile is somewhat larger than that of free-jet, 
the coefficient of contraction being correspondingly higher—i.e., Ce = (0.792) 
= 0.63.” At this point one would think that since the discharge coefficient de- 
pends upon the sum of contraction coefficients, the discharge coefficient also 
would be larger substantiating the experimental data given by Mr. McPherson 
and his associates for submerged discharge case. But Mr. Curtis continues 
to say that “Since the discharge coefficients of the two jets have already been 
seen to be essentially the same, (experimentally), one is forced to conclude 
that the rate of flow through the orifice is essentially independent of the tan- 
gential forces exerted upon the jet by the fluid into which it emerges.” This 
writer is forced to conclude that the departures of the characteristics of sub- 
merged flow case from free-discharge case requires further study and a bet- 


ter explanation. Unfortunately, the authors’ explanations can not be advocated 
as a panacea. 
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APPENDIX 


The integral given by equation (9) can be decomposed into two component 
integrals: 


at 


t dt a 
Vi )(t-a) (tor) 


Let us consider the transformation of Legendre:(12) 


2 
a + 


2 
c,sin + a, 


Ao, bo, Cg and dy can be so chosen that the above integrals simplify and re- 
duce to known integrals. It should be noted that these integrals, depending on 
the values of t, require separate treatment for each interval under con- 
sideration. 


I—In the intervals DE, and FEg, i.e., for t >d and t < -f, one has 


2 

def a(1+f 
(142) 


_4/(t-a) (142) 
sin y 0 <$.C¥2, 
(t-1)(def) 


2(d-1)(d+f)(14f) sin g,.cos g, 


[(aer)sin?g,- (142)] 


If Eqs. (15) and (17) are replaced in Eq. (13), one obtains 


t= (14) 
= 
and 
= (17) 
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2 
aM (def )sin“g,- a(1+f) ag, 
(19d) (det (142) 
aia a 
(18) 
(1¢f) (144) V1-ks sin“$, 


Integral (18) can be separated into a first and third kind elliptic integral: 


Q-= 1 a) F a d+f (19) 
{ (Py rks) + (a- -1) * 


7 9(9,n,k), is the third kind elliptic integral of Legendre, and may be written 
as 


Ths, on,k) 


(20) 


The parameter n-is taken as real, and may take any value, positive or nega- 
tive. For the intervals under consideration, 


<1 
Legendre’s integral is connected with Jacobi’s by the relation, 


= TL,(u, , & +1K’ (21) 


in which uy is written for F(Dj,k;). sn, cn, dn, are elliptic functions and K’ 
is the complete elliptic integral of the first kind for the complemantary modu- 
lus. €y is given by the relation 


an“ ) 


a+f asf 


In order to make the notations simpler, sn(€ 1,k1), cn(€4,k), and dn(é 4,k1) 
will be denoted by sn €j, cn€j, and dné€j, respectively. 
The solution of Jacobi’s integral is given by 


| 
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Ty(u, 2(E, + U, CnE, dné, H(é,-u) (22) 
2 H(€,+u) 


in which H(€ 1 * uy) is Jacobi’s Theta-function,(13) and 


(142) (tea 


Hence Legendre’s integral becomes: 


1 H(E,—u, 


cng, { H(é,¢u,) 


From Eqs. (19) and (23), we have 


H(€,- 
2 


which is the solution sought for the intervals under consideration. The con- 
stants M and N are determined from the comparison of 2-and t-planes for the 
corresponding values and were found to be one and ia resnectively. 

Il—For the interval HB, i.e., It! <1, 
If t is taken as 


-2f *) (25) 


2 (14f) 


and the foregoing manipulations repeated, one obtains 


2 
Qe= + (2-1) 
1+f 


V(14a)(142) 


in which 


=- which varies between zero and «2 

2 1+f 2 
Writing u, = F(g,,k,) and n, = - k2sn*(e ,,k,) = - <2 
2 2 2 a 1+f 


one has 


2 =| an) 


{ U5Z(E5,k5)+ 


e( EntBy ) 


2( f+a) 
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and therefore, 


in which 


g sin”* (1+f)(14t) 
2 


2( te?) 


a-1)(f-1 


sné = cn dn 


6 (€ tu) also is Jacobi’s Theta-function and connected to H(€ tu) by the relation 


@(z+ik’) = iBH(z) 


WK’ «iz 
4k 2K 


K and K’ are complete elliptic integrals of the first kind connected by the fol- 
lowing relationships: 


Il—For the interval BD, i.e., 1< t<d, 
If t is taken as 


(d-1 )sin”¢,, (14a) 


(29) 
(a-1)sin°g, - (144) 


one obtains 


N= 


Repeating the intermediate manipulations, one has 


where 
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in which 


x3 2( det) 
(142)(144) 


(a-1)(f£-1) ~ ‘2 2 


k = k «#k 


2 2 a-1 


= 
V (a-1)(ta1) 


IV—For the interval FH, i.e., -1 >t >-f 
If t is taken as 


a(f-1)sin°¥, - £(a41) 


(f-1)sin $4 + (de1) 


one obtains 


Bie { casa, (det) 


Viet) (144) 


in which 
f-1 2 2 


» sin@ = sn(é,,k4) 
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@(u,-1€, ) 
4 
an( €, 5%) @(u,s4€,) 


From Eqs. (33) and (34), one has 


Q= 2M ~(aet) sk, )-cn(E, ) 


(35) 


@(u,-L€ 


in which 


e 
u, = F( k,) = x? = —2(dst) 
* (142)(144) 


4 
(d+1) (£41) 


= 
fed 


and 


Intermediate steps in the solutions of the above integrals have been omitted 
for the sake of brevity. However, the reader will find an excellent account of 
the theory of Jacobian elliptic functions in reference (14), and in the other 
references stated previously. 

In order to obtain z from Eq. (11), dw should be written as, 


- +8 


The rest of the integration could only be performed by numerical methods. 
The best and well-known method for this purpose is King’s(15) A.G.M. scale 
method. 

If, for a fixed value of 8, numerical values are assigned to a, f, and d 
satisfying the boundary conditions in the transformed planes, corresponding 
Simultaneous values may be obtained for and anda. It should be 
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noted that this process was a tedious trial-and-error solution. For this rea- 
son, the plot in Fig. (4) was prepared to make possible the intermediate de- 
termination of an approximate result through the practical range of values. 
For large values ofa@, say greater than 15°, and for a given 8, Cc2 values 
are very nearly equal. Thus only three values of C-2 for a given B were 
computed and a curve drawn through them. Similarly for a fixed @ and fora 
wide range of values of 8, Cc1 values are very nearly equal. Therefore only 
four values of C, were sufficient to draw the appropriate curve. 

The writer believes that a more accurate numerical computation could be 
carried out. However, this has been precluded by the limited time available 
for the discussion. Besides, the above results are deemed to be sufficiently 
accurate for practical purposes. 


G. DUGAN JOHNSON, ! M. ASCE.—The authors do not concern themselves 
with the characteristics of, or problems associated with, butterfly valves for 
“wide open or fully closed” service in penstocks and at hydraulic turbine in- 
lets, but limit their discussion “to the flow of water through control-type 
valves for near-closed to near-open positions.” They say that “for most ap- 
plications flow characteristics are a function of geometry only” and attempt 
to demonstrate “the effects of geometric variables and present a basis for 
the analytical determination of flow coefficients” from which they indicate a 
method for predicting the inception of cavitation under regulating conditions. 

The authors state that “flow coefficients for the near-open positions are 
influenced by the relative blade thickness and the blade shape,” “viscous ef- 
fects are noted witha, = 0° to 20°” and “attainment of a high degree of test 
accuracy within this range is difficult because of the relatively large dis- 
charge rates and small head losses involved.” They also note that discrepan- 
cies in various published data for the near-closed positions “can be explained 
difference in closure angle.” 

For the intermediate valve disc positions, they developed a two- 
dimensional plastic flow section (8" wide by 1" thick) to observe and mea- 
sure the shapes of the two separate jets issuing from the orifices formed be- 
tween the side walls and the centerline profiles of the discs modeled. They 
reasoned that “since the jet area is a function of the variable (y), measure- 
ments of the jet thickness in a two-dimensional flow model, representing a 
narrow section of the center part of the blade, should provide a basis to un- 
derstand and predict prototype performance.” The writer agrees that such 
studies should be of qualitative assistance in understanding the characteris- 
tics of flow patterns at the centerline, but cannot accept the quantitative re- 
sults as being even approximately applicable to prototype performance. 

It is a well-known fact that the discharge and contraction coefficients of an 
orifice of any shape are critically affected by the edge condition, which must 
be accurately reproduced in the model in order to predict prototype per- 
formance. The authors recognized this when they added the dotted line to the 
right hand orifice in Figure 1, but they minimized it. unduly when they as- 
sumed that centerline flow patterns were even remotely similar to flow lines 
around any circular disc nearer its trunnions. 

Consequently, the writer feels that any attempt to predict prototype dis- 
charge performance for valves in pipelines on the basis of such two- 
dimensional model tests would only be adequate where extremely rough 
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approximations were acceptable, and that prediction of cavitation inception 
. .. . Which depends upon the maximum (rather than the average) velocity 
squared on the basis of such tests would be meaningless. 


HENRY VOLTMANN.!—The authors are to be commended for the very in- 
teresting, complete and valuable paper they have presented. Definite knowl- 
edge of butterfly valve performance becomes of vastly increasing importance 
in the hydraulic field and while many reports have been available of past 
tests, the data, as referred to by the authors, have not always been specific 
as to the test conditions. 

It may be interesting to add to the data presented the following information 
concerning a 6"' butterfly valve. 

Independent water flow tests were run on this butterfly valve in a “continu- 
ous pipe” installation by Mr. J. J. Colleville under the direction of Mr. M. B. 
McPherson at Lehigh University in 1953. The data from these tests agree 


very favorably with the data in Figure 7 and Table IV of the authors’ paper 
(nomenclature is identical): 


Rangeof Np No. of Test Runs 
5 
x 


x 


(Closure) - - 


The valve was lens-shaped with at eax /D of 0.188. There was no indica- 
tion of any’variation in the flow coefficient due to viscous effects. The C 
values obtained for individual test runs varied a maximum of 5% from the 
average values given in the above table, except at 70°, where only small rates 
of flow were attainable and the errors were consequently higher. 


The following piezometer station—distances were employed in the tests 
(four taps/station): 


Distance, in pipe diameters, from valve body 
Upstream 2D, 4D, 10D, 20D and 30D 


Downstream 1D, 2D, 3D, 4D, 5D, 7D, 9D, 11D, 13D and 15D. 


1. W. S. Rockwell Co., Fairfield, Conn. 


20° 1.00 2.5 - | 
30° 0.613 2.6 - 5.6 m 10° 6 7 
40° 0.38 2.5 - 5.6 x 10° 6 
45° 0.31 2.1 - 4.7 x 10° 5 - 
50° 0.23 0.5 -4.0 x 10° 11 
60° 0.13 1.3 - 2.7 x 10° 7 
70° 0.040 0.5 - 1.0 x 10° 5 ; 
80° 
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A straight piezometric head line was obtained for all tests with the up- 
stream stations, to and including the 2D station. The piezometric head line 
was straight in the downstream pipe, starting at the 5D station for blade 
angles of 10° to 40°, and starting at the 4D station for 45° to 70°. The valve 
loss, AH, was obtained by graphically extending the upstream and down- 
stream piezometric head lines to the valve and measuring the difference. As 


with the authors’ test data, the downstream pipe was filled with water in all 
test runs. 
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RONALD E. NECE.! J.M. ASCE.—The authors have presented the results 
of an interesting and comprehensive series of tests performed upon butterfly 
valves, and the information presented adds to the knowledge of this type of 
valve. The authors have also performed a service in presenting the results 
of some other investigators in a systematic form. Other data were undoubted- 
ly considered, but not presented, in this paper. Historically, the results of 
Weisbach, originally given in 1845 and tabulated in (D-1),2 were the earliest 
given for butterfly valves. 

Of the many types of valves used to regulate flows for different applica- 
tions in continuous pipelines, the butterfly seems to lend itself most readily 
to a theoretical analysis which will predict, with reasonable accuracy, the 
non-recoverable head loss through the valve. This discussion shall be con- 
fined to this one area of the investigations presented by the authors. 

The head loss caused by the presence of a constriction in a continuous pipe 
which causes a contraction of the flow, followed in turn by the rapid expansion 
of the stream to the full conduit area downstream from the valve, may be 
analyzed by the combination of energy and momentum considerations. 

This head loss may be expressed, using the authors’ notation, 


Ah 
29 


which is the familiar expression for the Borda-Carnot loss, where only the 
flow expansion downstream from the constriction is considered as being a 
source of loss. 


This expression may be rewritten as 
(2) 


. Asst. Prof. of Hydraulics, Dept. of Civ. and San. Eng., Massachusetts Inst. 
of Technology, Cambridge, Mass. 

. Numerals in parentheses refer to corresponding items in the list of 
References in this discussion. 
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where: C. = coefficient of contraction 


a = transverse projection of the restricted area of the valve. 

A = cross-sectional area of the valve body downstream from the 
constriction, in many cases identical with the area of the pipe. 

V = average flow velocity in the pipe. 


Equation (2) may be written in the common form 


2 
(3) 


where Ky = dimensionless head loss coefficient. For discharges generally 
associated with butterfly valves the Reynolds numbers are sufficiently high 

so that the flow geometry is essentially constant with varying discharges for 

a given valve configuration, with the result that K; is also a constant. Any 
theoretical expression for Ky thus depends upon the value of Ce; associated 
with the area ratio, a/A. 

The following simple “model,” shown in Fig. 1, may be considered as one 
method of use in predicting head loss characteristics. The valve is con- 
sidered to be in a two-dimensional duct of height B; the disk, which is pivoted 
at the duct center line, is of zero thickness, has the same width as the duct, 
and affects full closure of the valve at an angle a@ of 90°. The flow past the 
disk is then divided into two jets of widths Cc, b and Ccgab. 

For various angles a@ and slot height ratios 2b/B von Mises has calculated 
the values of Cc, for irrotational flow. These are listed in (D-2). From these 
calculated results the average value of C,, for the valve geometry may be ob- 
tained by interpolation for the assumed condition of equal pressures, and 
velocities, at the vena contracta of each jet. The distribution of the flow past 
the two ends of the disk is also determined as a part of this calculation. Thus, 
for the “model,” Ky; may be expressed as a function of @ alone. 

Limitations to the validity of the “model” are immediately obvious. The 
disks of actual valves possess a finite thickness of both blade and hub; thus, 
when the valve is fully open, there is still a reduction in flow area through the 
valve, with the consequent head loss. Actual disks are usually tapered toward 
the tips, and the tips sometimes additionally rounded; both of these conditions 
have the effect of tending to decrease the contraction of the streams on either 
side of the disk. For circular pipes and disks, the area “a” becomes a 
cresent-shaped opening, and the flow conditions are no longer two-dimension- 
al. Also, the disk angle at closure may be less than 90°; the 4"' and 6"" CDC 
valves reported by the authors had closure angles of 75° and 77.5°, respec- 
tively. The “model” also assumes no stops, recesses, or projections on the 
walls, and assumes that the tips of the disk, when at the full closure position, 
contact but are not recessed into, the valve walls. In addition, the two- 
dimensional measurements listed in the authors’ Table I indicate an unequal 
distribution of discharges around the two tips of the disk. 

Head losses predicted from the “model” would thus not be expected to fur- 
nish very good agreement with actual data at either extreme of disk angle. 
When the valve is fully open, the thickness of the disk and/or disk hub become 
predominate in the formation of the flow geometry; at, or near, full closure 
the calculated area ratio will obviously depart from the actual for valves hav- 
ing maximum @ values of iess than 90°. 
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In Table I are given the calculated values of Ky, for the “model” over the 
range of disk angles of 0°-90°. Also given in the table are the values of Ky 
obtained by the authors for the 4"' and 6"' CDC valves, values based upon test 
data by Dickey and Coplen and by Netsch and Schulz (authors’ references (10) 
and (9), respectively), as well as the older values given by Weisbach. The 

values of the discharge coefficient, Ca, as given by the authors, may be con- 
verted to those for Ky by equating comparable expressions for the head loss. 
This yields the relationship, 


(4) 


1.24 
Ki Ge 


The values of Ky for the authors’ references (9) and (10) were recalculated 
from the original papers. The values given are not precise, as they were de- 
termined from small-scale test data curves. 

As one of the series of check and control valves designed for specialized 
applications tested recently at the Hydrodynamics Laboratory of the 
Massachusetts Institute of Technology, some characteristics of the butterfly 
swing check valve “H” shown in Fig. 2 were determined. This valve differed 
from the others considered in this discussion in that there was a smooth 
transition from an 8" diameter pipe upstream from the valve to an 8"' square 
flow passage within the valve proper, followed again by a smooth transition to 
the circular pipe. The rectangular disk, allowed to rotate freely, was eccen- 
trically hinged as shown. For steady flows through the valve in the normal, 
or forward, direction, the angles @ assumed by the disk over the range of dis- 
charges of 1200 to 2000 GPM were 15° and 15°30’ for two slightly different 
shaped disks; the average values of Ky were 0.82 and 0.66, respectively. For 
discharges between 500 and 1200 GPM, the average value of Ky for the two 
disks were 0.65 and 0.72, respectively, indicating changes in the internal flow 
geometry with Reynolds number, these changes in this case possibly being 
caused more by the transition section effects upon the flow than by flow 
changes due to the valve disk. A description of the test facilities and results 
obtained for some of the other types of valves investigated are obtained in 
(D-3). 

Some of the data listed in Table I are plotted on Fig. 3. The solid line is 
designated as the “simple model theory.” The test points for swing check 
valve “H” are also plotted. It might be mentioned that data for some com- 
mercial valves of 10"' and 24"' diameter with circular disks also fall within 
the spread of the curves plotted. As expected, the differences between the 
“model” values and the test data are greatest for the largest and smallest 
values of @. 

The known physical size and shape of the valve and its disk for a particular 
valve could be used to determine more precisely the actual flow area ratios 
and to give better flow deflection angles in choosing values of C,. It is likely 
that by these more refined calculations the actual Ky, values for valves having 
a relatively simple interior could be estimated more precisely, dependent 
primarily, of course, upon satisfactory values of Cc. The simple “model” 
presented here is not to be considered a rigorous approach to the problem, 
and undoubtedly compensating factors in some cases cause its agreement with 
test data to be fortuitous. The complete characteristics of a given valve could 
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be determined only by a model test. However, the results clearly indicate 
that the butterfly valve is one type of control device where the application of 
the simple energy-momentum relationships will lead to a decent approxima- 
tion of head loss characteristics. The same approach may be used witha 
varying degree of success with other more common types of valves; i.e., gate, 
plug, etc. In smaller gate valves especially the head loss coefficient is found 
to vary with the size of the valve because of the relatively greater signifi- 
cance of the boundary configuration, gate recesses, and relatively thicker 
gate disk. Comparisons of various valves, using the area ratio a/A as the 
basis of comparison, are given in (D-3). 
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Discussion of 
“A STUDY OF BUCKET-TYPE ENERGY DISSIPATOR CHARACTERISTICS” 


by M. B. McPherson and M. L. Karr 
(Proc. Paper 1266) 


CORRECTIONS—As the result of a printer’s error, the illustrations that 
accompanied Proc. Paper 1266 were the illustrations that were used in Proc. 
Paper 1167. The illustrations for Proc. Paper 1266 are reproduced herewith: 
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CHARACTERISTICS OF BUCKET-TYPE ENERGY DISSIPATORS 
WITH ENTRANCE AND EXIT ANGLES OF 45° 


| 


0.1 


ik 
1348-60 HY 4 
| 
| 
08 | | , | 
hy | | 
| | | | 
YW) | 
=a 7 
0.4 05 0.6 0.7 08 0.9 Ho 
h 
FIGURE 


DISCUSSION 


SURGE HEIGHTS FOR BUCKET-TYPE ENERGY DISSIPATORS 
WITH ENTRANCE AND EXIT ANGLES OF 45° 
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SURGE HEIGHTS FOR BUCKET-TYPE ENERGY DISSIPATORS 
WITH EXIT ANGLE OF 45° AND ENTRANCE 
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